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ABSTRACT 


' A.  correlogram  is  the  two-dimensional 
autocoarrelati.cn  of  the  image  of  an  object 
i 3. lumi abated  with  non-coherent  radiation.  A 
laser  correlogram  is  obtained  from  the  power 
spectrum  of  the  irradiance  pattern  scattered 
from  fc be  object  when  illuminated  with  suf- 
fi.cien.-tly  coherent  radiation.  The  resolu- 
tion of  this  signature  is  dictated  by  the 
size  of  the  receiving  aperture,  with  rela- 
tively minor  degradation  by  atmospheric  tur- 
bulence. This  report  collects  in  one  place 
information  which  has  been  available  up  to 
now  only  in  conference  proceedings  or  in 
limited-circulation  Research  Notes  of  the 
Riverside  Research  Institute.  The  subjects 
treated  analytically  include:  a model  for 
laser  ToacRscattering,  studies  on  the  influ- 
ence of  atmospheric  turbulence  on  the  laser 
correLogram,  statistical  convergence  proper- 
ties o f the  laser  correlocram  signature,  qua- 
litative experimental  laboratory  results,  and 
the  outline  of  a design  for  a ruby-laser  ex- 
periment using  space  objects. 
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INTRODUCTION 

Ordinary  ground-based  imaging  of  space  objects  is  limited 
in  resolution  by  atmospheric  turbulence  to  a "seeing"  angle  of 
approximately  one  arc-second  (5  p,rad)  under  routine  clear-air 
observation  conditions.  One  method  of  circumventing  this  "see- 
ing" limitation,  laser  correlography.  can  yield  object  signa- 
tures with  "diffraction~limited  resolution".  This  Technical 
Report  is  a compendium  of  six  papers  which  go  into  varying  le- 
vels of  detail  on  several  aspects  of  laser  correlography.  Por- 
tions of  these  papers  have  appeared  previously  as  RRI  Research 
Notes1”5  or  in  conference  proceedings6*  7. 

The  "system  concept"  for  obtaining  ground-based,  diffrac- 
tion-limited signatures  of  space  objects  is  as  follows  (Paper 
No.  1) i The  signature  is  the  object  correlogram.  or  two-dimen- 
sional autocorrelation  of  the  object  radiance  in  the  absence 
of  turbulence.  This  signature  is  obtained  by  laser  correJog- 
raphv.  the  measurement  of  the  power  spectrum  of  the  irradiance 
pattern  scattered  from  an  object  when  it  is  illuminated  with 
sufficiently  coherent  radiation,  i.  e.,  when  it  lies  well  with- 
in the  coherence  volume.  Detection  of  the  backscattered  inten- 
sity across  the  receiver  aperture  is  employed:  a well-corrected 
receiving  telescope  is  not  needed;  when  fine  resolution  is  de- 
sired, a large  array  of  "photon  buckets"  can  be  used.  In  the 
special  case  in  which  a single  strong  glint  point  is  present  on 
or  near  an  otherwise  diffuse  target,  laser  correlography  is 
identical  with  the  well-known  on-axis  lensless  Pourier-trans- 
form  holographic  imaging  technique,6”10  In  the  absence  of  a 
glint  point,  a correlogram  is  obtained,  as  illustrated  in  Fig. 

4 of  Paper  No,  i.1*11*12  Systems  capable  of  yielding  useful 
signatures  in  practice  are  described,  and  energy  requirements 
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on  the  laser  illuminator  needed  for  space  objects  are  derived 
in  this  paper. 

Verification  of  the  basic  ideas  involved  in  laser  corre- 
lography  has  been  accomplished  at  RRI  by  a combination  of 
theoretical  work  and  laboratory  demonstrations.  Paper  No.  2 
includes  a model  for  laser  backscattering,  a condensed  account 
of  the  influence  of  atmospheric  turbulence  on  the  laser  corre- 
logram,  and  some  early,  qualitative  experimental  results.6 
These  results,  as  well  as  related  demonstration  experiments  per- 
formed by  others, e* 9* 11-14  are  consistent  with  our  calculations, 

A fuller  account  of  the  reason  why  laser  correlography  is 
insensitive  to  the  effects  of  atmospheric  turbulence  is  con- 
tained in  Papers  No.  3 and  4,  The  first  of  these  is  based  on 
the  Huygens-Fresnel  principle  and  a simplified  model  of  a rough 
target,  and  considers  an  idealized  laser  illuminator,  an  exo- 
atmospheric  space  object,  and  a ground-based  receiver,  with 
isoplanatism  assumed.  The  effects  of  using  a ground-based  il- 
luminator are  takeri  into  account  in  Paper  No.  4. 7 Realistic 
estimates  of  the  residual  effects  of  atmospheric  turbulence, 
related  to  scintillation,  are  given  as  a function  of  zenith 
angle.  For  slowly  rotating  space  objects  of  angular  subtense 
less  than  10  iirad,  it  is  concluded  that  the  system  concept  is 
valid  with  C02  laser  illumination  even  at  large  zenith  distan- 
ces if  a detector  array  is  employed,  and  with  a visible  laser 
of  long  coherence  length  for  the  object  very  close  to  zenith. 
Verification  of  the  concept  over  near-horizontal  paths  is  dif- 
ficult,7"’9 

The  statistical  convergence  properties  of  the  laser  cor- 
relogram  signature  are  examined  in  Paper  No.  5*  where  it  is 
shown  that  twenty  or  so  independent  realizations  of  the  irradi- 
ance  distribution  suffice  to  obtain  adequate  signal-to-noise 
ratio,5  This  result  has  been  verified  in  recent  demonstration 
experiments  at  RRI,15  (A  supplementary  Technical  Report  is  in 
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preparation15  which  contains  the  results  of  recent  laboratory 
demonstrations  and  computer  processing  of  laser  correlograms 
of  diverse  objects,  including  three-dimensional  "glinty"  ones.) 

Based  on  the  component  requirements  (suggested  by  Paper 
No0  1),  the  statistical  convergence  properties  (developed  in 
Paper  No.  5) t and  the  effects  of  atmospheric  turbulence  (Papers 
No.  3 and  4),  a feasibility  demonstration  experiment  was  de- 
signed. This  experiment  (Paper  No0  6)  would  employ  a ruby 
laser  transmitter  with  a coherence  length  of  10  meters,  opera- 
ting at  a repetition  rate  of  20  pulses  per  minute,  with  an 
output  of  10  joules  per  100-nanosecond  pulse.  A very  narrow 
(5  arc-sec)  beamwidth  is  required.  The  receiver  would  employ 
a SIT  or  EBS  tube  preceded  by  an  image  intensifier,  with  the 
S-25  photocathode  located  in  the  plane  of  the  (re-imaged)  en- 
trance pupil  of  the  (l.5-m  diameter)  telescope.  A full-fledged 
system  capable  of  imaging  high-altitude  space  objects  would 
employ  a C02  laser  and  a large  receiver  array,  as  discussed  in 
Paper  No.  1,  and  in  more  detail  in  Ref.  10 
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Paper  No.  1 

LASER  CORRBLOGRAPHY  SYSTEMS  FOR  OBTAINING 
SIGNATURES  OF  SPACE  OBJECTS 


(This  paper  is  an  edited  version  of  unclassified  portions 
of  Research  Note  N-5/I74- 3-33*  21  August  1972. 1 The  power 
calculations  relate  to  ensemble-averaged  power  within  a "single 
speckle".  For  other  (statistical)  considerations,  see  Paper 
No.  5.) 
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LASER  CORRELOGR&PHY  SYSTEMS  FOR  OBTAINING 
SIGNATURES  OF  SPACL  OBJECTS 


I.  Correlographv — An  Overview 

Correlography  is  a quasi -imaging  process  whose  output  is 
equal  to  the  autocorrelation  function  of  the  target  image; 
this  technique  may  provide  length  measurement  with  resolution 
limited  by  detection  aperture  in  addition  to  feature  and  target 
recognition.  A correlogram  system  detects  the  intensity  dis- 
tribution of  the  backscatter  ( self-interference)  pattern  of  the 
target  (Fig0  i)  and  operates  on  it  in  the  manner  shown  in  Fig.  2. 

The  average  intensity  of  the  backscatter  pattern  is  first 
subtracted  from  the  total  signal,  then  the  power  spectrum  of 
the  residue  (speckle  component)  is  computed,  followed  by  a 
two-dimensional  Fourier  transform;  it  has  been  shown  that  the 
resulting  two-dimensional  distribution  is  proportional  to  the 
image  autocorrelation  function;11  one  may  associate  with  the 
autocorrelation  function  a resolution  which  is  limited  by  the 
aperture  over  which  the  backscatter  intensity  was  detected. 

Such  correlogram  systems  are,  to  a large  extent,  immune  to  the 
atmosphere  because  they  operate  on  backscattered  intensity  and 
are  therefore  phase-insensitive.4  The  similarity  of  correlog- 
raphy to  holography  may  be  appreciated  by  noting  that  if  the 
image  consists  of  a large  diffuse  area  and  a bright  point  (such 
as  that  due  to  a glint)  then  the  correlogram  is  identical  to 
the  on-axis  hologram8  °. 

Using  the  setup  of  Fig.  3>  correlograms  were  formed  for 
simple  shapes  in  the  laboratory,  yielding  the  results  shown  in 
Fig.  4.  The  top  row  shows  the  target  images  that  were  used  to 
generate  the  correlograms;  the  second  row  contains  the  correlo- 
grams of  the  images;  the  third  row  contains  traces  formed  by  a 
pinhole  detector  scanned  through  the  principal  axes  of  the  cor- 
re3ograms  of  the  rectangle  and  T shapes.  Several  features  and 
properties  of  correlograms  become  apparent  upon  inspection  of 
this  figure.  First,  it  is  seen  that  the  correlogram  offers 
information  along  two  dimensions  which  are  both  transverse  to 
the  optical  line-of-sight  and  that  these  dimensions  may  yield 
accurate  measurements  of  lengths  of  various  target  features. 
Secondly,  all  correlograms  have  point  symmetry,  imposing  some 
ambiguity  on  the  patterns;  this  may  make  unique  target  identi- 
fication difficult  in  some  cases.  Finally,  the  correlograms 
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In  Fig.  1,  the  schematic  diagram  of  a ground-based  laser 
system  which  serves  as  a model  for  the  discussions  of  this 
paper*  the  relative  scales  of  the  laser  illuminator  and  the 
elements  in  the  receiver  array  are  distorted.  Also  note  that 
the  laser  illuminator  would  typically  be  placed  somewhere 
within  the  receiver  array.  It  is  implicit  throughout  the 
discussion  that  meteorological  conditions  permit  line-of-sight 
optical  viewing  of  the  space  object;  no  heavy  cloud  cover*  fog* 
haze*  or  precipitation  is  presumed  to  lie  within  the  path0 
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FIG.  2 CORRELOGRAPHY 


IF  I (x)  = IMAGE, 

IMAGE  ACF  = j*l(x)l(x+u)dx 


IMAGE  AUTO- 
CORRELATION 
FUNCTION  (ACF) 

CORRELOGRAM 
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are  observed  to  preserve  angle  relationships  among  paddles, 
aerials,  and  other  protruding  linear  structures. 

In  view  of  the  potential  value  of  these  signatures,  RRI 
has  made  a study  of  means  and  rational  procedures  for  infor- 
mation retrieval  from  correlograms  of  interesting  targets. 5* 6* 15 
This  study  includes  mensuration,  feature  recognition,  target 
recognition,  and  image  recovery.  In  addition,  laser  correlo- 
grams have  been  formed  in  the  laboratory  for  interesting  tar- 
gets for  the  purpose  of  revealing  the  nature  of  the  degrada- 
tions caused  by  such  effects  as  the  three-dimensional  nature 
and  multiple  glint  content  of  typical  targets.15 

II0  Atmospheric  Limitations  (See  also  Papers  No,  2,  3*  and  4.) 

It  is  known  that  propagation  of  a laser  illumination  beam 
through  the  atmosphere  to  a space  object  is  influenced  by  the 
turbulence  structure  of  the  atmosphere;  in  particular,  it  ap- 
pears that  the  beam  will  break  up  and  this  may  cause  a highly 
structured  illumination  to  fall  on  distant  targets,  resulting 
in  severe  fading  of  the  returns. 7 In  addition,  the  intensity 
distribution  of  the  backscatter  pattern  is  known  to  be  per- 
turbed by  a random  apodization  due  to  the  atmospheric  structure; 
this  random  apodization  may  limit  the  ultimate  resolution  of 
laser  backscatter  systems,  providing  an  upper  limit  on  useful 
aperture  sizes  for  receiver  arrays. 4* 6,7 

It  is  recommended  that  the  influence  of  the  atmosphere  on 
target  illumination  and  signature  resolution  be  explored  in 
order  to  place  ultimate  limits  on  performance  of  laser  back- 
scatter systems.2 

III.  System  Design  and  Implementation 


Possessing  a general  understanding  of  the  signatures  offer- 
ed by  laser  correlography  systems,  one  can  consider  the  diffi- 
culties envisioned  in  implementing  practical  operating  systems. 
This  section  discusses  some  of  the  principles  and  specifications 
for  the  selection  of  the  illumination  laser,  backscatter  detec- 
tion system,  target  tracking  system,  and  receiver  configuration. 

The  discussion  is  presented  with  reference  to  the  general 
system  block  diagram  shown  in  Pig.  5*  The  major  elements  in 
the  system  are  the  system  control  computer,  tracker,  illumina- 
tor, array  receiver,  and  data  recording  system.  The  system 
control  computer  processes  coarse  acquisition  data  and  fine 
tracking  data  in  order  to  point  the  illumination  laser  and  the 
array  elements;  it  also  provides  synchronization  of  illuminator 
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pulses,  receiver  range  gates,  and  the  data  recording  system. 

A separate  target  tracker  is  provided  to  convert  the  coarse 
acquisition  data  (which  may  have  been  obtained  from  an  LWIR 
system)  into  fine  tracks  with  accuracy  on  the  order  of  10  jirad. 
The  tracker  consists  of  a tracking  laser,  duplexed  transmit/ 
receive  optics,  an  angle  and  range  tracker  receiver,  and  a 
mechanical  pointing  system.  After  an  adequately  precise  track 
is  established,  the  illumination  subsystem  directs  ■‘'he  illumi- 
nation laser  beam  through  the  transmit  optics  toward  the  target; 
a precision  pointing  system  is  required  to  keep  the  illumina- 
tion beam  on  the  target.  Backscattered  signals  are  detected 
by  an  array  of  receiver  telescopes;  each  array  element  consists 
of  a precise  pointing  system,  receiver  optics,  and  a range- 
gated  receiver  for  detecting  the  backscattered  energy.  The 
data  recording  system  stores  all  pertinent  backscatter  signals 
and  system  diagnostic  data. 

A.  WaveJ.emth  Utilization 

In  this  section,  preferred  operational  wavelengths  are  es- 
tablished, general  design  concepts  are  proposed,  and  the  speci- 
fication of  the  illumination  and  array  receiver  subsystems  are 
discussed. 

Figure  6 compares  the  "diffraction-limited  resolution" 
(taken  as  the  product  of  laser  wavelength  and  target  range,  di- 
vided by  array  aperture  length)  for  ruby  and  carbon  dioxide 
lasers  as  a function  of  receiver  aperture  size,  with  target 
range  as  a parameter.  Note  that  with  ruby  laser  illumination, 
a 1.5-meter  aperture  can  yield  resolution  of  20  cm  for  400~km 
target  range  and  50  cm  on  1,000-km  targets;  carbon  dioxide  sys- 
tems with  like  resolution  would  require  apertures  of  20  m and 
50  m,  respectively.  For  10,000  km  target  range,  a C02  laser 
correlography  system  requires  an  aperture  of  400  m in  order  to 
obtain  resolution  of  40  cm  on  a target;  proportionally  larger 
apertures  are  required  for  synchronous  altitudes.  It  is  con- 
cluded that  it  appears  more  practical  to  employ  visible  laser 
wavelengths  in  order  to  reduce  the  large  aperture  requirements 
associated  with  carbon  dioxide  systems. 

However,  Fig.  7 compares  ruby  and  carbon  dioxide  trans- 
mitter energy  requirements  as  a function  of  the  number  of  pho- 
tons received  per  spatial  resolution  cell  in  the  backscatter 
pattern  for  altitudes  up  to  several  thousand  kilometers,  (The 
derivation  of  these  curves  rs  given  in  the  Appendix.)  It  is 
easily  observed  that  visible  lasers  require  much  more  trans- 
mitter energy  to  yield  the  same  number  of  photons  in  the  recei- 
ver. (it  is  shown  in  the  Appendix  that,  at  low  altitudes,  the 
required  laser  energy  is  proportional  to  the  cube  of  the  optical 
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FIG.  6 RESOLUTION  VS  APERTURE  SIZE  FOR  TWO  WVELENGTHS 
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FIG. 7 ENERGY  REQUIREMENTS  FOR  LOW-ALTITUDE 

LASER  CORRELOGRAPHY 
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frequency.)  For  instance,  a target  at  1,000  km  altitude  requires 
2 mJ  from  a C02  laser  and  10  J from  a ruby  laser,  respectively, 
to  return  100  photons  per  resolution  cell.  On  the  basis  of  this 
energy  comparison,  it  is  concluded  that  visible  lasers  are  not 
likely  to  have  adequate  energy  to  provide  useful  signal-to-noise 
ratios  for  target  ranges  in  excess  of  roughly  1,000  km  and  that 
C02  lasers  are  most  suitable  for  the  synchronous  altitude  role. 

B.  Visible  Wavelength  System  Parameters  (See  also  Paper 

No.  6.2) 

The  objectives  of  a visible  wavelength  measurements  program 
would  be  to:  l)  evaluate  the  quality  of  the  correlography  and 
holography;  2)  confirm  single-  and  multiple-glint  analyses; 

3)  obtain  10-20  cm  data  on  known  and  unknown  targets;  4)  test 
computation  schemes;  and  5)  combine  these  data  with  information 
obtained  from  other  kinds  of  sensors. 

Table  I compares  typical  current  and  modified  specifications 
of  a ruby  laser  for  making  such  measurements;  note  that  state- 
of-the-art  pulse  energy  is  adequate,  but  that  both  the  pulse 
length  and  the  coherence  length  must  be  extended  to  make  them 
longer  than  the  physical  lengths  of  targets  ejected  to  be 
observed. 

Fig,  8 is  a schematic  diagram  of  the  receiver  configuration 
needed  to  record  the  backscatter  patterns  for  computerized  analy- 
sis of  the  signatures.  The  field  lens  images  the  intensity  dis- 
tribution in  the  aperture  of  the  telescope  onto  a photoelectric 
image  tube;  the  intensity  pattern  on  the  image  tube  (seen  through 
a band-oass  optical  filter)  is  then  read  into  a recording  system 
where  off-line  data  processing  is  performed  to  extract  target 
information  ,> 

C.  Carbon  ID.oxide  System  Parameters 

Figure  9 shows  the  transmitter  energy  needed  to  obtain  vari- 
ous signal-to-noise  ratios  in  each  receiver  element  for  target 
altitudes  from  2,000  to  4,000  km;  a 1-m  receiver  aperture  is 
assumed  along  with  10-dB  path  loss,  10-p,rad  illumination  beam- 
width,  and  a target  reflectivity-area-product  equal  to  1 m2. 
(These  curves  follow  from  the  energy  considerations  presented 
in  the  Appendix. ) Figure  10  relates  to  the  maximum  allowable 
duration  of  the  laser  illumination  pulse  such  that  the  back- 
scatter pattern  moves  less  than  one-half  a speckle  size  during 
the  illumination  time.  It  is  seen  that  for  target  angle  rates 
of  up  to  10  mrad/sec,  pulse  lengths  of  50-100  [xsec  are  adequate 
for  target  lengths  of  up  to  10  m transverse  dimension. 

Because  of  the  large  aperture  sizes  that  are  needed,  the 
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TABLE  L RUBY  LASER  SPECIFICATIONS 


SPECIFICATION 

STATE-OF-THE-ART 

MODIFIED 

PULSE  ENERGY 

10  J 

10  J 

PULSE  LENGTH 

15  NSEC 

*100  NSEC 

COHERENCE  LENGTH 

1-2  cm 

10  M 
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FIG.  10  MAXIMUM  TARGET  ANGLE  RATE  ALLOWABLE  FOR 
VARIOUS  LASER  PULSE  LENGTHS 
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requirement  for  array  receivers  (as  opposed  to  contiguous  aper- 
tures) is  self-evident  for  carbon  dioxide  systems.  The  number 
of  elements  required  in  a filled  array  (an  array  with  a uniform 
grid  of  equally  spaced  telescopes,  where  the  spacing  is  equal 
to  one-half  the  minimum  expected  speckle  size)  may  be  related 
to  the  number  of  resolution  cells  desired  on  the  target. 

Pig.  11  is  a plot  of  the  required  number  of  array  elements  as 
a function  of  linear  resolution  on  the  targets,  for  maximum 
target  dimensions  of  3 and  10  m.  It  is  clear  that  10-cm  reso- 
lution on  a 10-m  target  will  require  10,000  elements,  and  that 
this  is  unacceptable.  The  conclusion  is  that  receiver  array 
thinning  is  essential;  in  addition,  the  thinning  must  eliminate, 
for  reasons  of  cost  and  practicality,  over  99^  of  the  elements 
in  order  to  reduce  the  numbers  to  a maximum  of  100-200  elements. 

A sketch  of  a receiver  subsystem  that  would  be  suitable  for 
energy  detection  in  each  array  element  is  shown  in  Fig.  12. 
Doppler  tracking  is  performed  by  an  electronically  tracking 
local  oscillator  after  an  optical  heterodyne  has  been  performed 
on  the  photo-mixer  element.  A square-law  detector  detects  the 
energy  in  the  IP  passband  and,  after  low-pass  filtering,  the 
signal  is  recorded.  An  alternative  receiver  scheme  might  make 
use  of  a tracking  optical  local  oscillator,  although  this  tech- 
nology is  less  adequately  developed.  This  receiver  is  applic- 
able to  on-axis  holography  and  correlography;  for  off-axis 
holography,  using  polarization-sensitive  reception,  the  incoming 
optical  signal  would  first  be  passed  through  a polarizing  beam 
splitter  and  then  be  transmitted  to  two  superheterodyne  receivers 
operating  in  the  same  way  as  the  one  shown  up  to  the  IF  output. 

At  that  point,  the  output  of  the  two  IP  amplifiers  would  be 
added  and  passed  through  a square-law  detector;  the  low-pass 
filtered  output  would  contain  the  required  interference  between 
the  two  bedims.  The  polarization-sensitive  receiver  would  re- 
quire coherent  optical  and  electronic  loca1  oscillators  in  each 
channel  to  operate  successfully. 
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FIG.  II  NUMBER  OF  ARRAY  ELEMENTS  FOR  EQUAL 
SPACING 
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APPENDIX 

LASER  ILLUMINATOR  ENERGY  REQUIREMENTS 


If  a pulsed  laser  beam  of  angular  width  Rj  , energy  W.  and 
wavelength  X is  directed  on  a target  of  projected  area  a,1 
diffuse  reflectivity  p and  range  R,  then  the  average  number  of 
photons  backscattered  per  unit  area  in  the  vicinity  of  the 
transmitter  is  given  by 

n = W-j  X p ox  a , (1) 

1:  2ir  he  R4 


where  h is  Planck's  constant,  c is  the  speed  of  light  and  xa  is 
the  round-trip  atmospheric  transmission  factor. 

In  order  to  preserve  the  structure  of  the  backscatter 
pattern,  the  aperture  area  of  the  receiver*  must  be  smaller  than 
one-half  to  one-fourth  the  smallest  possible  area  of  a speckle 
on  the  ground,  the  latter  quantity  being  given  by  X2/fit,  where 
is  the  expected  maximum  solid  angle  subtended  by  a target; 
for  10-m-diameter  targets,  this  limits  aperture  diameters  to 
less  than  1 m for  ranges  up  to  about  30,000  km  for  ruby  and 
2,000  km  for  carbon  dioxide  systems.  If  it  is  assumed  (somewhat 
arbitrarily)  that  maximum  aperture  diameters  for  each  element 
of  an  array  receiver  will  be  limited  to  1 m by  issues  of  cost 
and/or  practicality  then  it  may  be  seen  that  two  regimes  of 
photon  integration  area  nay  be  defined;  in  the  low-altitude 
regime  the  area  over  which  photons  are  integrated  is  limited 
by  speckle  size,  and  in  the  high-altitude  regime  this  area  is 
limited  by  receiver  element  aperture  area.  In  the  low-altitude 
regime,  receiver  photon  integration  area  is  given  by 


Ar  = KX2/fit  = KX2R2/At  (LOW  ALTITUDE)  (2) 


where  At  is  the  minimum  circular  area  required  to  inscribe  the 
maximum  target  dimension  and  where  K~l/2  to  1/4;  for  high- 
altitude, 

Ar  = CONSTANT.  (HIGH  ALTITUDE)  (3) 

*Or,  more  precisely,  the  area  over  which  photons  may  be 
integrated,  and  properly  called  a detection  resolution  cell. 
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The  transmitter  energy  required  to  produce  a given  number 
of  photons  Nr  at  the  photodetector  for  receiver  area  Ar  may  be 
computed  from  £q.  1 and  2 or  3 for  both  altitude  regimes  after 
rioting  that  the  optics  transmission  factor  t0  may  be  combined 
with  xa  to  yield  a total  transmission  factor  t.  Thus, 

W-  = 2ir  hc  niA4-R2Nr  (LOW  ALTITUDE)  (4) 

1 K A 3 p ot 

or 

w.  _ 2tt  hc  ft-jR^Nr  (HIGH  ALTITUDE)  (5) 

1 AporAj. 

Equations  4 and  5 were  used  in  connection  with  Figs.  7 and  9, 
respectively.  (N  = nr  Ar>) 
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Paper  No.  2 

INFLUENCE  OF  THE  ATMOSPHERE-  ON  LASER-PRODUCED  SPECKLE  PATTERNS 

(This  is  an  unaltered  reproduction  of  a paper  prepared  for 
the  1974  International  Optical  Computing  Conference,  April  9-lij 
I974,  Ztirich,  Switzerland. 6 The  paper  includes  a model  for 
laser  backscattering,  a condensed  account  of  the  influence  of 
atmospheric  turbulence  on  the  laser  correlogram — a fuller 
account  being  contained  in  Paper  No,  3 — and  some  qualitative 
experimental  laboratory  results.  More  recent  results  are  con- 
tained in  a forthcoming  Technical  Report.15) 


{ 
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INFLUENCE  OF  THE  ATMOSPHERE  ON  LASER  PRODUCED  SPECKLE  PATTERNS* 


by 

M.  King,  M.  Greenebaum  and  M.  Elbaum 
Riverside  Research  institute 
New  York,  N.Y.  10023 


We  analyze  the  situation  of  radiation  leaving  a coherently  illuminated  rough  object  and 
passing  through  the  turbulent  atmosphere  to  a plane  of  observation  where  the  speckle  pattern 
is  measured.  Employing  a form  of  the  Huygens -Fresnel  principle  which  applies  to  i spatially 
inhomogeneous  random  medium,  we  find  that  the  effects  of  the  atmospheric  turbulence  may  be 
represented  as  a random  apodization  of  the  speckle  pattern  which  would  be  obtained  in  .the 
non-turbulent  case,  provided  that  the  entire  target  lies  within  a single  isoplanatic  region. 
Using  experimental  dataon  stellar  scintillation  to  estimate  the  degree  of  random  apodization 
for  a ground-based  receiver  observing  a near-zenith  exo-atmospheric  target,  we  estimate  the 
strength  of  apodization  at  selected  wavelengths.  This  analysis  shows  that  one  can  expect 
minimal  depth  of  the  modulation  at  10  pm;  however,  patterns  formed  in  the  visible  region  of 
the  spectrum  are  inclined  to  suffer  substantial  apodization.  Nevertheless,  depending  on  the 
scale  size  of  "atmospheric  speckles"  relative  to  the  range  of  speckle  sizes  produced  by  the 
target,  speckle  pattern  parameter  estimation  in  the  plane  of  observation  may  still  be 
possible.  This  possibility  is  demonstrated  by  means  of  a laboratory  simulation  experiment. 
Thus,  at  least  in  principle, so  long  as  the  normalized  autocovariance  function  of  the  random 
apodization  is  known,  both  the  power  spectrum  and  autocorrelation  of  the  speckle  pattern  of 
the  object  may  be  determined  directly  from  the  statistics  of  the  combined  irradiance 
distribution;  these  statistics,  in  turn,  yield  the  correlogram  (autocorrelation)  and  power 
spectrum  of  the  object. 


1.  Introduction 

It  is  well  known  that  the  speckle  pattern 
generated  by  a coherently  illuminated  object 
carries  specific  information  about  the  object; 
i.e.,  the  power  spectrum  of  the  speckle 
pattern  is  equal  to  the  autocorrelation  func- 
tion of  the  object  illumination.  1*2  In  this 
paper,  we  evaluate  how  well  this  irradiace 
power  spectrum  can  be  estimated  when  it  is 
subject  to  degradation  by  atmospheric  turbu- 
lence . 

First,  we  shall  derive  a formula  for  the 
power  spectrum  of  the  speckle  patterns  formed 
from  the  coherent  radiation  scattered  by  an 
object  (two-dimensional,  for  simplicity),  and 
propagated  through  the  atmosphere.  This 
formula  will  serve  as  a basic  tool  in 
evaluating  the  strength  and  the  nature  of  the 
atmospheric  perturbation.  The  validity  of 
some  of  the  conclusions  derived  analytically 
is  demonstrated  by  experiments  described  at 
the  end  of  the  paper. 

2.  Formal  Results 

A coherently  illuminated  object  is 
described  here  by  a field  distribution 
EQ  ($;t).  By  using  the  tangent  plane  approx- 
imation, this  field  can  be  written  as 

Ec(S,t)  — t (?)  l0%©,t)exp[i  k s(?)]f(l) 


where  I (?,t) is  the  illumination,  t is  the 
reflection  coefficient,  s (?)is  a random 
variable  describing  the  surface  depth  variation 
about  the  mean  plane,  andk  is  a wave  number. 
For  rough  objects, it  is  customary  x'/to  assume 
that  s(p)is  Gauss ian-distributed,  with  mean 
zero,  and  standard  deviation  a not  smaller 
than  the  wavelength  X of  the  radiation  used 
to  illuminate  the  target.  As  a result  of 
a s s X tne  autocorrelation  function  of  the 
object  Jrield  distribution  can  bo  well-approxi- 
mated by  a delta  function.  Consequently,  the 
fields  from  different  points  at  the  object  are 
independent,  and  the  distribution  of  this 
field  in  space  can  be  regarded  as  a stochastic 
process  with  white  noise-like  spectrum. 

The  field  at  a Fraunhofer  observation 
plane  of  the  object  is  computed  by  using  3 
the  extended  Huygens-Fresnel  principle, 

E (p‘ , t)  E(p.,t)  G (&,?’),  (2) 

0 i=l  1 

where  N is  the  number  of  independent  "cells" 
intv.  which  the  field  at  the  object  was  de- 
composed, ana  G is  the  Green’s  function 

G (^,3’)-  |exp[ikh  (?.,?’)  + Y (?.,2*)]  (3) 

= Gq  (?..$•)  exp  [*  ©i,**)],  (3a) 

where  z is  the  distance  between  the  object 
plane  and  the  receiver  plane;  the  phase  term, 
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T l (3^3')  +10*  (3^3’ > 

is  a complex  phase  perturbation  which  de-  , 
scribes  the  effects  of  the  inhomogeneous 
medium,  0®  and  4 representing  the  phase  and 
the  logarithm  of  the  amplitude  perturbations, 
respectively,  due  to  the  random  inhomogeneir 
ties  of  the  atmospheric  index, 

Gq  (3^3’)“  I exp  [i  It  h^,?’)] 

is  the  value  of  the  Green's  function  in  tha 
absence  of  the  turbulent  medium,  and  h(p,,p  } 
is  the  geometric  distance  between  the  points 
and  p'. 

The  "instantaneous"  irradiance  of  the 
observation  point  $ , !(?*)#  i®  obtained  by 
time-averaging  the  square  of  Re  S E (p ' , t)  f 
over  times  small  compared  with  the  time  it 
takes  for  any  to  change: 

I(3')  = H Re  jE(3',t)  E*  <3’,t)  } (5) 

- h E0*<fy 

Go(3i,3,)  GQ*  (3j *3' (3^3') 
+Y*(3i,3,)]|  <5*> 

= >5  exp  2-1  (3r?')Reg^E0(3i)E0(?*j) 

G0  (3^3')  Gq*  (3j.P*)}  . (5b) 

where  the  Eq.  (5a)  is  obtained  by  inserting 
Eqs.  (3)  and  (4)  into  Eq.  (5). and  Eq.  (5b) 
results  from  the  assumption  that  the  object 
lies  well  within  the  isoplanatic  region,  i.e., 
that  Y(?.,3')  = Y (PVP')  with  i^j  and  i,j  « 
1,2. . . ,N.  BeginningJwith  Eq.  (5a),  the 
time-dependence  of  the  fields  is  suppressed 
to  simplify  the  notation, 


In  the  absence  of  the  random  medium, 

Y = 0 , 1(3')  reduces  to  1’0(3‘)  » the 
object  speckle  pattern: 


I0(3')  = h Re  }E  (3*)  Es* 

where  Es(3‘)  =^1E0^i^  G0(3i»3')*  (6) 


X (3‘),  the  expected  value  of  I (3')»  in  the 
sense  of  an  ensemble  average  over  many  real- 
izations of  the  surface  roughness,  can  be 
shown  to  be  equal  to 


I (3’)  = 


i=l 


(7) 


The  term,  IA(3').  is  referred  to  as  the 
atmospheric  speckle  noise.  It  is  the 
irradiance  distribution  which  would  be 
present  at  the  receiver  point  3'  when  a 
coherent  point  source  of  strength  T (3')  is 
present  at  some  object  point  pQ: 

Ia(3!)  » X (3')  exp  [2  4(3'  , PQ)]  , (8) 

It  has  been  shown4  that  the  expected 
value  of  I,  taken  as  an  ensemble  average 
over  many  A realizations  of  the  atmospheric 
turbulence  is  equal  T (3')>  Using  Eqs.  (8) 
and  /£)) , we  rewrite  the  expression  for  the 
instantaneous  irradiance  distribution  I(3')> 
in  terms  of  the  expected  value  of  I(p'), 

X (3'),  randomly  modulated  by  both  the  object 


and  the  atmospheric  epeckles,  i.e.. 


i&’Kxfp') 


1(3’)  -I 

\ , Xo(3')  - i(p')  1 

L K3‘)  J . 


(9) 


It  ia  easy  to  aee  that  the  normalized  autg- 
covarience  function  of  I(3‘)»  C(3',p'  + Ap’), 
can  be  written  as  a sum  of  three  terms  : 
a target-related  term  which  is  the  normalized 
autocovarience  function  of  I (3 ' ) , an  atmo- 
sphere turbulence -related  term,  which  is  the 
iiormalizedautocovariar.ee  function  of  IA(3’)> 
and  their  product. 

c(3',?'  + Ap' ) = =^Jj  <p-  ± AP'I-^V 

I 


= CA  + C0  + CAC0  ' 

tfher<1  0a  . 


(10) 


_ <Xp  (3‘)  +Ap')>  t 


I2 


and  1(3')  = 1(3'  + Ap')  = I has  been  assumed. 


The  spectrum  of  I (3) , S,  is  computed 
by  taking  the  Fourier  transform  of  the 
normalized  covariance  function  (10) 


S = Ss  + S + S * S , 
A O A o 


(ID 


where  SA  and  SQ  are  the  Fourier  transforms 
of  Ca  and  Cq  respectively,  and  * denotes  . 
the  convolution  operation. 


3.  Discussion  of  Formal  Results 

• We  want  to  know  how  well  the  spectrum 
S of  the  speckle  pattern  randomly  modulated 
by  the  atmospheric  noise  approximates  the 
object-related  power  spectrum,  SQ. 

Several  conclusions  can  be  drawn  from 
a slightly  modified  form  of  Eq.  (10) : 

C = CQ  + CA  (CQ  +1)  . (10a) 

We  observe  first  that  the  function  C approxi- 
mates the  function  of  interest,  C0»  if  at 
least  one  of  the  following  conditions  is 
satisfied: 

(a)  C.  < < c s 1 for  all  values  of  the  lag, 

a3'.  0 

(b)  « h where  and  L are  measures  of 
the  Width  of  the  correlation  functions 

<IA  (^')  IA  (3’  + A3')>  and  < I0(3,)I0(3’+A3'» 

respectively. 

The  inequality  (a)  requires,  in  partic- 
ular, that  the  normalized  standard  deviation 
of  turbulence  generated"atmQspheric  speckle 
noise", CA  (3'/3')  = a(  / T2  be  much  smaller 
than  the  normalized  standard  deviation  of  the 
object  related 'bignal,"  Cs  (p',3')=:  Og  / I “1 
Inequality  (a)  therefore  signifies  that  the 
depth  of  modulation  of  the  atmospheric 
speckles  considered  as  a"noise  "is  much 
smaller  than  that  of  the  "signal". 
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On  the  other  hand, inequality  (b)  is 
satisfied  when  the  half-width, La,  of  the 
correlation  function  of  the  "noise", 

< *A  (?')  *a  ($'  + A£')  > , is  much  smaller 
than  the  half-width,  L0  , of  the  correlation 
function  of  the  " signal", <ID (p ' )I0{p'4Ap' )>  t 
A3  one  can  see  from  Eq.  (10a)  the  function 
C (j$',jS'  + Aj$') closely  approximates 
C0 ($',?•  + Ap’jfor  all  values  of  the  lag 
&p’  larger  than  La,  when  (b)  is  satisfied.  In 
the  “noise "-related  second  term  of  Eq.  (10a), 
Ca  imposes  a "lag  window”  on  the  "signal"- 
related  term  CQ.  The  narrower  the 
correlation  function  C»  with  respect  to  the 
correlation  Cc,  the  wider  the  spectrum  of 
"noise"  Sa  will  be  with  respect  to  the 
spectrum  of  "signal"  Sc.  If  the  favorable 
situation  of  condition  (b)  exists,  one  may 
scan  the  noisy  irradiance  distribution  with 
an  aperture  small  enough  to  resolve  the 
speckle-pattern  "signal"  and  sufficiently 
large  to  average  out  the  turbulence -generated 
"noise",  (i.e.,  integrate  over  the  structure 
of  characteristic  size  La) . 

We  may  now  turn  to  the  available 
experimental  data  in  order  to  find  out  under 
what  conditions  the  atmosphere  is  sufficiently 
cooperative  for  (a)  to  hold.  It  is  known 
from  stellar  scintillation  data  that,  for 
vertical  propagation  in  the  visible  spectrum, 
the  normalized-varjance  of  the  turbulence- 
related  term  is  approximately  "0.2. 

The  situation  at  10.6  u,at  least  for  vertical 
paths  at  night  and  scaling  the  visihle 
scintillation  according  to  the  X-7'6  law, 
is  expected  to  be  even  more  encouraging, 
i.e.,  **  0.008.  However,  for  near- 

horizontal propagation  under  "saturation-of- . 
scintillation"  conditions,  e.g.,  propagation 
over  paths  in  excess  of  '=>_1  km  in  the  visible 
or  " 10  km  at  10.6  um,  a!l  /T2  «,  10  and  the 
condition  (a)  is  not  lA  satisfied.4-6 

As  to  the  second  condition,  (b) , we 
note  that  the  width, L,  of  the  autocorrelation 
function  of  the  "signal"  is  proportional  to 
Xz/D,  where  D is  the  "aperture  size"  o"  the 
object,  e.g.,  the  diameter  of  a uniformly 
reflecting,  uniformly  illuminated  object.  On 
the  other  hand,  the  "noise"  correlation 
length,  LA,  is  approximately  (XL)*5,  where  L 
is  on  the  order  of  the  path  length  through 
uniform  turbulence  (e.g.,  horizontal 
propagation)  or  the  height  to  the  tropopause 
for  vertical  propagation.5' %ased  on 
measurements  of  ^shadow  bands"in  stellar 
astronomy,  (X  L)^  can  be  expected  to  vary 
from  3 or  4 to  7 or  8 cm  in  the  visible.  6 


* The  third  condition,  L « Lft,  is  likely 
to  be  of  practical  importance  unless 

condition  (a)  is  also  satisfied.  For  details, 
see  Ref.  4. 


j.  Condition  (b) , which  requires  that 
(XD^  < Xz/D,  is  satisfied  for  objects  with 
transverse  dimensions  of  meters,  at  ranges 
10  km  or  greater  in  the  visible. 

In  summary,  the  theoretical  prediction 
supported  by  experimental  data  shows  that 
there  are  situations  where  it  will  be  possible 
to  obtain  reliable  information  about  the 
object  from  the  power  spectrum  of  the  object 
speckle  pattern  randomly  apodized  by  the 
atmosphere.  These  situations,  given  above 
as  conditions  (a)  and  (b) , are  more  easily 
satisfied  in  the  infrared  than  in  the  visible 
spectrum. 


4.  Experimental  Results 


To  simulate  a diffusely  reflecting  object, 
an  apertured  ground-glass  diffuser  was  used 
to  phase-modulate  a collimated  laser  beam 
(X  = 0.5  ^m)  in  transmission.  _ A silicon-diode 
array  vidicon  (effective  aperture  0.13  mm  x 
0.01  mm)  located  1.8  ra  from  the  object,  was 
capable  of  resolving  the  minimum  speckle  size 
in  the  speckle  pattern  scattered  by  the  object. 
Directly  in  front  of  the  diode  array  («  3 mm 
away)  wa3  placed  a photographic  transparency 
containing  a simulated  "atmospheric  speckle 
pattern".  (This  transparency  was  made  by 
recording  the  speckle  pattern  produced  by  a 
2 x h nun  2 diffuse  target  at  a distance  of 
1.8  meters.  After  4 min  of  development  in 
DK-50  at  68  ° F the  mean  amplitude  transmission 
of  the  Ektapan  recording  film  was  approximately 
50%.  Consequently  the  irradiance  fluctuations 
cause  a severe  amplitude  modulation  to  be 
recorded  which  subsequently  appears  on  the 
field  passing  through  the  transparency.  The 
film  backing  and  the  emulsion  thickness 
variations  tend  to  produce  some  phase  modu- 
lation which  must  be  taken  into  account  in 
the  detailed  interpretation.) 


Vhe  output  of  the  vidicon  provides  an 
ar.alog  spectrum  analyzer  with  an  irradiance 
distribution  (e.g..  Fig.  1 (b) ) which  is  used 
to  obtain  the  power  spectrum  of  the  sample. 
Individual  power  spectrum  estimates  were 
obtained  by  sampling  in  one  direction  at  0.01- 
mm  intervals  over  7.5  mm  of  the  speckle  pattern 
and  averaging  over  0.13  mm  in  the  orthogonal 
direction. 


To  obtain  a sufficient  number  of  random 
samples,  in  order  to  estimate  the  power 
spectrum  to  an  empirically  sufficient  degree  of 
accuracy,  31  different  realizations  of  the 
object  phase  structure  were  generated  by 
exposing  non-overlapping  areas  of  the  ground 
glass  to  the  collimated  beam.  The  arithmetic 
mean  of  these  individual  spectral  realizations 
(displayed  as  Sp  and  SD',  respectively  in 
Figs.  (1)  and  (2)  ) estimates  the  power  spectrum 
of  the  object  speckle  pattern.  Thi3  can  be 
compared  with  the  estimator  of  the  power 
spectrum  { S and  S'  in  Figs.  (4)  and  (5), 
respectively)  of  the  object  speckle  pattern 
randomly  apodized  by  the  photographic  trans- 
parency. (The  spectra  in  Figs.  (4)  and  (5) 
are  obtained  as  an  arithmetic  mean  of  31 
individual  spectrum  estimates,  each  obtained 
for  a particular  realization  of  the  object 
phase  structure  and  of  the  film  transparency 
perturbation.) 
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Of  the  two  objects,  one  (Fig.  1 (a)) 
had  a wider  spectrum,  and  the  other  (Fig.  2 
(a))  a narrower  spectrum,  than  that  of 
the  simulated  atmosphere  (Fig.  3).  Upon 
comparirig  the  spectra  displayed  in  Figs.  1 
and  4 and  those  in  Figs.  2 and  5,  we  observe 
that  the  general  features  of  the  spectra  are 
preserved  in  the  presence  of  the  perturbation 

i.e.,  neither  size  nor  shape  is  seriously 
affected. 

The  authors  wish  to  thank  L.  Schlom 
and  P.  Wong  for  carrying  out  the  experiments. 


* This  work  was  supported  by  ARPA. 
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Fig.  3:  (a) "Atmospheric"  speckle  pattern, 
(b)  "S  ",  Spectrum  of  (a)  vs. 
spatial  frequency,  f^. 

T-l/306-3-11 


Fid.  5:  Spectrum  of  speckles  from  object 
in  Fig.  2(a)  as  modulated  by 
"atmospheric"  speckles  in  Fig. 3 (a) 


Paper  No. 


LASER-PRODUCED  SPECKLE  PATTERNS  AS  SEEN  THROUGH 


A RANDOMLY  INHOMOGENEOUS  ATMOSPHERE 


(The  following  is  a reproduction  of  the  latest  revision 
of-RRI  Technical  Paper  PTP-IOO(r),  21  December  1973.  Previous 
versions  of  this  paper  appeared  under  the  same  title  as  Techni- 
cal Paper  PTP-100,  25  June  1973,  and  as  Research  Note  N-8/174- 
3-33,  28  November  1972. 4 ) 

This  paper  is  a first-principles  account  of  the  reason 
why  laser  correlography  is  insensitive  to  atmospheric  phase 
fluctuations,  and  includes  some  estimates  of  the  residual  ef- 
fects of  turbulence  (related  to  scintillation).  A ground-based 
receiver,  an  exo atmospheric  space  object  with  a rough  surface, 
isoplanatism,  and  an  idealized  (space-based)  laser  illuminator 
are  assumed.  See  Paper  No.  ^ for  the  effect  of  relaxing  some 
of  these  constraints. 
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LASER-PRODUCED  SPECKLE  PATTERNS  AS  SEEN  THROUGH 
A RANDOMLY  INHOMOGENEOUS  ATMOSPHERE 

Michael  Greenebaum  and  Marvin  King 
Riverside  Research  Institute,  80  West  End  Ave.,  New  York,  N.  Y.  10023 

The  spatial  autocovariance  function  of  the  irradiance 
distribution  received  in  the  backscatter  direction 
from  a target  illuminated  uniformly  by  a laser  is  con- 
sidered in  the  presence  of  an  intervening  turbulent 
atmosphere.  When  the  target  is  perfectly  diffuse  ana 
lies  entirely  within  a single  isoplanatic  region  of 
the  receiver  aperture,  the  effect  of  the  atmospheric 
turbulence  is  to  produce  a random  apodization  of  the 
irradiance  distribution.  Prom  this  randomly  apodized 
irradiance  distribution,  one  may  recover  the  auto- 
covariance function  which  would  be  obtained  in  the 
absence  of  turbulence.  This  non-turbulent  autoco- 
variance function  is  proportional  to  the  far-field 
diffraction  pattern  corresponding  to  the  object  illu- 
mination function,  while  the  associated  power  spectral 
density  is  proportional  to  the  target  "correlogram" . 
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Introduction 

It  is  we31  known  that  the  backscattered  irradiance  from  a 
rough,  laser-illuminated  target  has  a speckled  appearance. 1-4 
The  second-order  statistics  of  this  irradiance  distribution 
convey  spatial  information  concerning  the  target  in  a manner 
similar  to  the  van  Cittert-Zernike  theorem:  the  spatial  auto- 
covariance function  of  the  random  irradiance  pattern  is  propor- 
tional to  the  far- field  diffraction  pattern  corresponding  to 
the  object  illumination  function;  the  associated  power  spectral 
density  is  proportional  to  the  "correlogram"  of  the  object 
illumination  function. 2-4  In  the  present  paper,  we  examine 
the  autocovariance  and  power  spectral  density  of  the  irradiance 
for  the  case  where  a turbulent  atmosphere  intervenes  between 
the  target  and  the  entrance  pupil  of  the  receiving  telescope. 

To  keep  the  underlying  interference  phenomena  in  the  fore- 
front, we  base  our  analysis  on  the  random-phase,  identical- 
strength,  point-scatterer  target  model  and  calculational  ap- 
proach of  Goldfischer.3  Using  a form  of  the  Huygens-Fresnel 
principle  which  applies  to  a spatially  inhomogeneous  random 
medium, 5-6  we  extend  this  analysis  to  the  case  where  the  radia- 
tion leaving  a coherently  illuminated  target  passes  through 
the  turbulent  atmosphere.  Provided  that  the  entire  target  lies 
within  a single  isoplanatic  region,  the  effect  of  the  atmospheric 
turbulence  is  shown  to  reduce  to  a random  apodization  of  the 
irradiance  distribution  which  would  be  obtained  in  the  non- 
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turbulent  case.  The  auto covariance  function  is  represented 
therefore  as  the  sum  of  a target-related  speckle  term,  a term 
arising  from  the  random  apodization  by  the  turbulence,  and 
their  product  (Eq.  (20)). 

When  experimental  data  on  stellar  scintillation  are  used 
to  estimate  the  degree  of  random  apodization  to  be  expected, 
we  conclude  that  the  turbulence  effect  is  small  in  favorable 
cases,  e.  g.,  a ground-based  receiver  observing  a near-zenith 
exoatmospheric  target  illuminated  at  10.6  [xm.  Then  the  second- 
order  statistics  retain  their  "van  Cittert-Zernike  theorem-like" 
significance.  In  general,  it  is  possible  to  extract  the  desired 
target-related  second-order  statistics  rom  the  randomly  apodized 
irradiance  distribution,  regardless  of  the  strength  of  the 
apodization,  so  long  as  the  normalized  autocovariance  function 
of  the  apodization  is  known  from  auxiliary  data  (e.  g.,  on 
appropriately  filtered  starlight). 

Although  the  details  are  not  carried  out  explicitly  here, 
the  methods  used  in  this  paper  are  capable  of  fairly  straight- 
forward generalization  to  estimate  the  effects  of  violating  the 
isoplanatic  patch  condition. 

The  coordinate  system  shown  in  Pig.  1 is  used  throughout. 
Points  p'  on  the  coherently  illuminated  target  have  coordxnates 
(u,v, O).  The  radiation  scattered  from  the  (exoatmospheric)  tar- 
get passes  through  the  turbulent  atmosphere  to  the  (ground-based) 
telescope  receiver  aperture  located  at  a distance  h.  The 
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irradiance  at  each  point  p in  the  entrance  aperture  is  recorded. 
The  points  p have  coordinates  (x,y,h).  For  simplicity,  a 
planar  target  is  assumed,  and  the  plane  of  the  receiver  aper- 
ture is  taken  as  parallel  to  the  plane  of  the  target. 


1 
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The  Generalized  Huvqens-Fresnel  Principle 

For  propagation  distances  large  compared  with  the  size  of 
a radiating  aperture,  it  has  been  shown5 * 8 that  the  Huygens- 
Fresnel  principle  may  be  extended  to  the  spatially  inhomogene- 
ous random  medium.  According  to  this  principle,  the  field  due 
to  an  arbitrary  complex  disturbance  over  the  radiating  aper- 
ture can  be  computed  by  superposition  of  the  fields  due  to 
sph  ?rical  wavelets  leaving  each  aperture  element.  The  propa- 
gation (f  an  arbitrary  wavefront  is  then  determined  by  the  way 
a spherical  wave  propagates  through  an  instantaneous  realiza- 
tion of  the  inhomogeneous  random  medium.  The  complex  field 
at  the  observation  point  p due  to  a point  source  at  p'  m the 
plane  of  the  transmitting  aperture  (which  is  identical  to  the 
field  at  p'  due  to  a point  source,  at  p)  is  given  by:' 

G(p‘,p)  = fl(p,p»)  = s"1  exp  Ciks  +^(p,p,)H  * U) 

where  s = j p - p’  | is  the  geometric  distance  between  the 
field  point  p and  the  source  point  p’, 
k = 2tt/X  = optical  wavenumber,  and 

yCPiP1)  = -t(p,p' ) + i /z5a(p,p»)  describes  the  effects 
of  the  inhomogeneous  medium  upon  the  propagation 
of  a spherical  wave.6 

Re  ij/  = l represents  the  logarithm  of  the  amplitude  perturbation 
due  to  the  random  inhomogeneities  of  the  atmospheric  refractive 
index,  and  Im^  = represents  the  phase  perturbation  arising 
from  the  same  inhomogeneities. 
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Calculation  of  Backscattered  Irradianee 

As  the  first  step  in  our  calculation  of  the  second-order 
statistics  of  the  irradianee  distribution,  we  consider  the  com- 
plex. field  scattered  from  the  target  when  the  propagation  takes 
place  in  a non-turbulent  medium.  The  target  is  modelled  as  a 
collection  of  diffuse  scatterers,  each  occupying  an  area  Au  Av 
« ?v2  over  which  the  reflected  irradianee  P(u,v)  = P(u,v,  0)  may 
be  considered  constant.7  Each  scatterer  contributes  a field  at 
the  observation  point  p at  time  t,  which  is  dependent  upon  the 
reflected  irradianee  P(u,v)  leaving  the  target  at  the  point 
p«  = (u,v,  0)  at  the  time  t*  = t - s/c,  and  upon  /rfs(u,v),  a 
random  phase  angle  associated  with  the  scatterer  at  p* . 

Now  we  let  a spatially  inhomogeneous  medium  intervene 

— & 

between  p*  and  p.  By  means  of  the  extended  Huygens-Fresnel 
principle,  Eq.  (l),  we  may  obtain  the  field  contribution  at 
time  t at  point  p = (x,y, h)  due  to  the  same  scatterer  after  the 
wave  has  propagated  through  the  turbulent  atmosphere,  i.  e., 

[tP(u,v)  Au  Av/rrS2H  1//a  Re  e*P  D-k(ct  - s)  + i /rf*(u,v)Il 

...(2) 

where  c - speed  of  light  in  the  medium, 

^(u,v)  = /rfs(u,v)  - i^(p,p’ ) = /rfS(u,v)  - i^(u,v,0;x,y,h), 

T = atmospheric  transmission  between  (u,v,0)  and  (x,y,h), 
P(u,v)  = reflected  irradianee  (or  "object  illumination 
function")3,  the  product  of  the  irradianee  incident  upon  the 
target  surface  and  p^(u,v),  the  (hemispherical)  diffuse  reflec- 
tance at  (u, v).  P(u,v)  is  the  "scene". 
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The  correlation  between  the  real  random  variables,  /$s(u,v),  and 
the  complex  random  variable,  ^(u, v, 0;x,y, h),  is  identically 
zero,  since  these  quantities  are  associated  with  physically 
distinct  and  independent  processes.  As  discussed  in  Appendices 
A and  B,  /5s(u,v)  is  assumed  to  be  Gaussian-distributed,  with 
zero  mean  and  standard  deviation  greater  than  tt,  8j  9 whereas 
^(u,v,  0;x,y,  h)  is  assumed  to  be  a complex  Gaussian  process  with 
effectively  uncorrelated  real  and  imaginary  parts,  an  assump- 
tion which  is  well-founded  experimentally. 10 

In  the  Fresnel  approximation,  the  contribution  of  the 
scatterer  at  (u,v, 0)  to  the  field  at  (x,y, h)  is  given  by: 

OrP(u,v)  Au  Av/wh2!]  x/a  exp  Q,(u, v;x,y)H 

1 • cos  C\(ct  + 2£U±XV)  + G(x,y)  + jrf‘s(u,v)  + j4a(u, v;x,y)H  , 

...(3) 

where  0(x,y)  = - ^(h  + ) > and 

^«s(u,v)  = j*S(u,v)  + *^(u2  + v2)  . 

The  weak  trend  in  the  phase,  /^,s(u, v),  produced  by  the 
quadratic  Fresnel  term  associated  with  the  scatterer  at  (u,v, 0) 
may  be  absorbed  without  significantly  affecting  the  random- 
ness.3 Thus,  ^,s(u,v)  is  also  a Gaussian-distributed  random 
phase  with  large  variance.8  The  z-dependence  of  jrfa(u, v;x,y) 
is  suppressed  in  the  notation;  f6a  gives  explicitly  the  phase 
perturbation  introduced  by  the  atmospheric  inhomogeneities 
which  are  present  along  a particular  path,  from  (u,v, 0) 
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to  (x,y,h).  Slriiilar  comments  hold  for  the  notation  given  the 
log-amplitude  perturbations,  l{  u,  v;  x,  y) . 

E(x,y)  =*  E(x,y,h),  the  total  field  at  (x,y,h)  in  the  re- 
ceiver aperture  (Fig.  l),  is  obtained  by  superimposing  the 
contributions  from  all  scatterers  in  the  (u,v, 0)  plane: 

B(*,y)  = V2  {=os(^ct  + 9(x,y) ) £ 

•exp  P-t(u, v;x,y)H  ccs  C X'U^y~)  + ^,s(u,v)  + /rfa(u, v;x,y)H 

- sin(^ct  + 0(x,y))  ^2J  [>(u,  v)!]1/2 
•exp  Q,(u,v;x,yO  sin  C ^-(XU^YV)  + ^’s(u,v)  + *$a(u,  v;x,yO  } 

We  now  average  the  square  of  E(x,y)  over  times  long  com- 
pared with  the  period  of  the  optical  radiation,  but  short  com- 
pared with  the  time  it  takes  for  /£a(u,  v;x,y),  <t(u,v;x,y),  or 
^'s(u,  v)  to  change.  The  result  is  the  "instantaneous"  irradi- 
ance  at  (x,y,h): 

H( x,y)  = { (£  &(u^v)^2  ex?  C-t(u,vjx,y)I] 

V 

•cos  C ^(XU^~)  + /^,s( u,v)  + ^a(u,v;x,y) 

+ ( 2 CP(u,v)lK2  exp  C-t(u,  v?x,y)II 

.sin  C ^*u±*V)  + ^|3(u,v)  + flu,V!x,y) 

•••(5) 

The  sum  of  the  squares  of  each  term  in  the  series  is  inde- 
pendent of  the  phase  angles  of  the  individual  scatterers,  as 


— liO— 
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well  as  of  the  turbulence-related  phase  angles.  This  sum,  which 
is  given  by 


5(x,y)  = 


£ P(u,v)  exp  E”2£(u, v;x,y)3 
u,  v 


T P pw 

Au+O*  lih*  J du  I dv  p(“»v)  exp  Dat(u,vjx,yO  « 
Aw0  -*  -*  • ...(6) 
is  the  irradiance  distribution  which  would  be  obtained  under 
noncoherent  illumination  of  the  target,  for  a single  realization 
of  the  random  apodisation  produced  by  the  atmosphere.  In  other 
words,  Eq.  (6)  may  be  considered  to  result  from  an  average  over 
an  ensemble  of  random  target  phases,  jt$s(u,v),  for  a single 
realization  of  the  turbulence. 

The  terms  remaining  in  the  series  expansion  of  Eq.  (5)  after 
Eq.  (6)  is  subtracted  will  be  denoted  by  the  symbol,  H(x,y), 

1.  e.,  H(x,y)  = H(x,y)  - H(x,y).  In  the  case  of  a uniformly 
illuminated  target  and  no  turbulence  in  the  intervening  medium, 
H(x,y)  represents  the  variation  of  the  speckle  pattern  about  the 
mean  over  the  plane  of  observation.3  More  generally,  in  the 
presence  of  turbulence,  we  interpret  H(x,y)  as  the  difference 
between  the  irradiance_measured  at  the  plane  of  the  receiver 
aperture  when  the  target  is  coherently  illuminated  and  the  ir- 
radiance which  would  be  measured  upon  scattering  from  a similar 
noncoherently  illuminated  target.  (The  same  instantaneous  reali- 
zation of  the  turbulent  atmosphere  and  of  the  rough  surface  is 
presumed  under  both  illumination  conditions.)  Expanding  Eq.  (5) 
and  deleting  all  terms  contributing  to  H(x,y),  we  obtain: 
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H(*,y)  - ^ 


E E &(u,v)  p(u  - v - 

u,v  <D,fi 


• exp  £t(u,v;x,y)  + l(u  - v - |j^;x,y)II 

•J"cos  Owx  + fjy  + ¥(u),ft;u,V7X,y)[3  - 6^0  •J  • (7) 


where 


- « - if . v'  ■ v - if . *• 

CD  = |g(u  - U»)  , n - |g(v  - V»)  , 
¥(co,a;u,v;x,y)  « /$,s(u,v)  - /rf*s(u 

+ /rfa(u,v;x,y)  - ^“(u  - 

_ i 1 when  a = b, 

6ab  *"  * 0 when  a ^ b. 


_ 2aM  v _ 2&£k\ 

2tr  ' V 2n 


v 


*hft  ..  . 

2T;A,Jf;  * 


and 


The  notation  used  is  similar  to  that  of  Ref.  3?  o>/2rr  and  Q,/ 2n 
are  spatial  frequencies  in  the  x-  and  y-direct.ions,  respectively. 
The  Kronecker  deltas  remove  the  terms  for  which  a>  = ft  = 0,  ones 
that  contribute  to  H(x,y).  Any  single  term  in  Eq.  (7)  describes 
the  interference  fringe  generated  in  the  plane  of  observation  by 
a pair  of  point  scatterers  on  the  diffuse  surface,  one  at  the 
point  (u,v,  0)  and  the  other  at  the  point  (u  - •jj—,  v - 0), 

as  seen  through  the  intervening  randomly  inhomogeneous  atmo- 
sphere. v?x,y)  is  a random  phase  angle,  so  that  one  has 

a superposition  of  a random  collection  of  high-contrast  inter- 
ference fringes,  so  long  as  the  observation  time  is  short  com- 
pared with  the  time  it  takes  for  <t(u,  v;x,y)  and  ^a(u,  v;x,y)  to 
change  appreciably. 11  Moreover,  it  is  clear  from  the  definition 
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that  Y(0, 0;u,v;x,y)  s 0,  for  any  point  (u,v)  on  the  target  sur 
face  and  any  point  (x,y)  in  the  telescope  entrance  aperture. 

Letting  angular  brackets  denote  the  ensemble  mean,  i.  e., 
the  average  over  all  possible  realizations  of  the  random  vari- 
ables, -t,  /rfa,  and  ^s,  we  show  in  Appendix  A that 

< H(x,y)  > = 0.  ( 

It  follows  that  the  ensemble  mean  irradiance, 

< H(x,y)  > = < H(x,y)  > 


= J du  J dv  P(u,v)  < exp  2<t(u,v;x,y)  > 

■ J du  J av  “ s°  » 


—00  —CO 


the  mean  value  of  the  irradiance  in  the  absence  of  turbulence, 
since  conservation  of  energy12  requires 

< exp  2-t(u,v;x,y)  > = 1.  ( 


T-1/306-3-H- 


RIVERSIDE  RESEARCH  INSTITUTE 


Spatial  Autocovariance  of  the  Irradiance 

From  the  expressions  for  the  irradiance  H(x,y),  Eqs.  (6) 

and  (7),  we  can  obtain  directly  the  spatial  autocorrelation 

function,  < H(x,y)  H(x+£,y+£)  >,  i.  e.,  the  ensemble  average  of 

• • 

the  product  of  the  irradiances  at  pairs  of  points  in  the  re- 
ceiver aperture.  If  we  assume  that  the  random  variable  K is 
effectively  uncorrelated  with  the  random  variables,  f6a  and 
(in  the  sense  required  in  Appendix  A),  then  in  the  notation3 

Aco  = Au  ; Av  , (ll) 


the  autocorrelation  of  the  total  irradiance  distribution13 
becomes: 


< H(x,y)  H(xfr$,y+C)  > = yg“£rg  ^ ^ p(u'> v’ )HX/2Au  Av 

J-Ott  n u^v  u.^v, 

• E E [>(u  - |^,v  - £h-%  p(u*  - AteL,v.  - ^-0  */»  Aa>  Afl 
cD»,n»  n n 


»<exp  [t(u,v;x,y)  + <t(u',v»  ;**•£, y+C) 

+ ^(u  - ~^,v  - ^?x,y)  + £(u»-  ^v'-  ^^-?3c*-|,S+C)H  > 


2tt 


2it 


2n 


.<(  cos  Qux  + Qy  + I (cd,£2;u,v  ;x,yD 

.cosQo'(xH)  + «'(y+C)  + !(<»', fl'ju'.v'ix+^y+CO)  • (12> 

It  is  shown  in  Appendix  B that  the  last  factor  in  the  typical 
term  of  the  above  series  expansion  reduces  to: 

6 6 . 6n,  + i Re  fexp  i(u£+fiC)  «<  exp  iA^a(x,y?£,C)  > 

(1)0  fio  (D'o  Si'o  C L _ 

•Ej6  6 6 6 + 6 "xivn  ^ TvhH  ^ _ 

u,u'  v,v’  a),a)’  u’,u-  v»,v-  -5^“  cd,-<d» 


2rr 
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where 

A>Sa(x»y = - ^a(u(v;x,y)  + ^a(u  - |^,v  - ^:x,y) 

* Au.vjxH.y+C)  - )la(u  - |jj“  v - (13b) 

is  a phase  angle  which  depends  upon  the  differences  in  the  atmo- 
spheric phase  perturbations  at  pairs  of  points. 

The  first  term  in  Eq.  (±3a)picks  out  the  "d-c"  part  (at 
zero  spatial  frequency,  co  = Cl  = co*  = = 0)  of  Eq.  (12).  Upon 

comparison  with  Eq.  (6)  and  the  change  of  variables,  Eq,  (ll), 
this  d-c  part  is  seen  to  be  equal  to  the  autocorrelation  of 
H(x,y),  i.  e.: 

< H(x,y)  H(x+£,y+C)  > = £ Au  Av  P(u,v) 

" u,  v 


• £ Au  Av  P(u',v»)  < exp  [>t(u, v;x,y)  + 2-t(u» , v>  ;xH,y+C)I]  > 


u’ , v' 


Au  *■»  0 W^h4  * 
Av->  0 


J du^j  dv  P(u,v)  exp  [^(u,  v;x,y)H 


i*00  r00  —I 

-ooj  du-®J  dv  p(u>v)  ex&  L.2't(u,v?^,4,y+OJ 


> 

...(14) 


The  remainder  of  Eq.  (13)  consists  of  two  terms  which  yield 
equal  contributions13  when  the  summations  over  u',v’,^S  and 
are  performed  in  Eq.  (12) . Together,  these  contributions  repre- 
sent the  autocorrelation  function  of  H(x,y),  i.  e. : 


< H(x,y)  E{y&i,y+0  > 

= < H(x,y)  H(x+£,y+£)  > - < H(x,y)  H(x+£,y+£)  > 
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r Re  j du  j dv  r*  f°dQ  exp  i(o>£  + fl£) 


lfcm 


—00  <—00  —00  —eo 


< exp  iA/*a(x,yj£,e)  > P(u,v)  P(u  - |jp,v  " 


2tt 


. < exp  &(u,v;x,y)  + *(u  - 

+ <t(u,v;x+4,y+£)  + ^(u  - ^,v  - > 

(The  first  equality  in  Eq.  (15)  follows  from  the  result, 

< H(x,y)  H(xt4,y+0  > = < H(x,y)  H(xt^Y+0  > = 0,  . 


...(15) 


which  one  can  show  readily.) 

In  the  absence  of  turbulence  (l  0,  A/rfa  + 0),  Eq.  (15) 
expresses  the  well-known  results  that  the  variance  of  the 
irradiance  is  equal  to  the  square  of  the  mean  irradiance, 13 
<(  H(x,y)  ys  = Ho2,  as  required  by  the  negative-exponential 
irradiance  statistics,2  and  furthermore,  that  the  covariance 
function  of  the  irradiance  distribution  is  proportional  to  the 
two-dimensional  Fourier  transform  of  the  target  autocorrelation 
function, 2 4 i.  e.,  the  far-field  diffraction  pattern  corres- 
ponding to  the  object  illumination  function.3 

When  turbulence  is  present,  we  consider  first  the  impor- 
tant special  case  that  the  entire  target  lies  within  a single 
isoplanatic  region. 14  Then  we  may  assume  that  for  any  point  p 
on  the  receiver  aperture,  the  complex  phase  perturbation, 

^(p,p»)  = 't(p,p‘)  + i ^a(p,p*)  is  the  same  for  all  points  p* 
on  the  target.  Hence,  we  may  take  A/z$a(x,y;£, C)  = 0 in  Eqs.  (13) 
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and  (15).  Furthermore,  although  £(u,v;x,y)  may  indeed  depend 
upon  the  (x,y)  coordinates  of  the  receiver  aperture,  we  may 
take  «t(u,v;x,y)  = -t(u  - ^~,v  - |~;x,y),  etc.,  in  Eqs.  (13), 

(l4)  and  (15),  outside  the  integrals.  Thus,  the  autocorrelation 
function  of  the  "noncoherent  irradiance",  H(x,y),  reduces  to: 

< H(x,y)  H(xt£,y+£)  > 

= < e*P  [>t(u,v;x,y)  + 2£(u,  v;x+£,y+£)]]  > Hq2,  (16) 

where  (u,v)  are  the  coordinates  of  a typical  point  on  the  target, 
and  H , the  mean  irradiance  in  the  absence  of  turbulence,  is 
given  by  Eq.  (9).  Similarly,  the  spatial  autocorrelation  func- 
tion  of  H(x,y)  - H(x,y)  - H(x,y)  becomes: 

< H(x,y)  H(xt^,y+C)  > 

= < exp  H2^(u,v;x,y)  + 2£(u,  v;xK;,y+£)I]  > CHs(x,y;x+|,y+C), 

...(17) 

where  cHs(x,y;3d-e,y+(;)  = < H(x,y)  S( , y+C )>non_turtoulent 

T 2-s  2 P°°  P°° 

* l^A^  Re  J ^ J exp  i(ai^0 
• J au_J  dv  p(u,v)  p(u " *5 v ■ Ir* 

is  the  spatial  autocovariance  of  the  irradiance  fluctuations  in 
the  absence  of  turbulence. 3) 13 

The  common  multiplicative  factor  in  Eqs.  (16)  and  (17) 
arises  from  the  random  apodization  of  the  irradiance  by  the 
atmospheric  turbulence,  the  same  phenomenon  which  produces 
scintillation. 10  A similar  scintillation-related  random  apodiza- 
tion (and  similar  independence  of  atmospheric  phase  perturbations) 
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has  been  noted  in  the  analysis  of  Fourier "transform  holography 
in  the  presence  of  a random  medium.  14“le  The  degree  of  random 
apodization  to  be  expected  in  typical  situations  in  the  visible 
and  infrared  regions  will  be  estimated  in  the  next  section. 

Finally,  if  the  entire  target  does  not  lie  within  a single 
isoplanatic  region,  the  ensemble  averages  in  Eqs.  (±4)  and  (15) 
must  be  evaluated  before  the  integration  is  performed.  Similar 
ensemble  averages  arise  in  estimating  the  atmospheric  modulation 
noise  in  an  optical  heterodyne  receiver17  and  in  estimating  the 
size  of  the  isoplanatic  region  in  Fourier-transform  holography. 15 
Although  we  omit  the  details  here,  it  is  clear  that  when  the 
target  extends  beyond  a single  isoplanatic  region,  one  does  not 
ejqject  to  obtain  the  neat  factoring  into  the  product  of  a target- 
related  part  and  a turbulence-related  part  which  one  finds  in 
Eqs.  (16)  and  (17)  for  the  autocorrelation  functions. 

Returning  to  the  isoplanatic  situation  for  which  Eqs.  (l6) 
and  (17)  apply,  let  us  evaluate  the  common  factor, 

Fa  s < exp  D*t(u,v;x,y)  + 2l{u, v;x+£,y+£)H  > 

= < exp  2£(u, v;x,y)  >.<  exp  2-t(u,  v;xt£,y+C)  > 

* exp  ^(u,v;x,y;x+?,y+0, 

= exp  4c^(u,v;x,y;jd-4,y+C),  (18) 

where  the  spatial  autocovariance  of  the  random  function  f(x,y) 
is  given  by  Cf(x,y;x+£,y+£)  = < Cf(x,y)  “ < f(*,y)  >3 

• |jE(x+£,yK)  - < f(x+£,y+0  >'J>. 
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The  second  equality  in  Eq.  (18)  follows  from  the  assumption  that 
the  log-amplitude,  -t(u,v;x,y),  is  ?.  Gaussian  random  variable. 10# 12 
C.  (u, v;x,y;jd-£,y+£)  is  the  log-amplitude  covariance  function 

v 

for  a point  source  at  (u,v).  The'  final  equality  follows- from  Eq.  (10). 
Alternatively,  F may  be  expressed  in  terms  of  the  spatial 

a 

autocovariance  function  of  the  irradiance  distribution  received 
in  the  (x,y)  plane  from  the  single  point  source  at  (u,v), 

CHa(x,y;xt£,y+C)  “ < H(x,y)  H(x*i,y+C)  >point 

source 

- < H(x,y)  > < H(x*£,y+C)  > 

= Hq2  < exp  2<t(x,y)  ><  exp  2<t(x+£,y+C)  > 

• J^exp  4c^(u,v;x,y;^,y+C)  - l] 

= Hq2  (Fa  - i).  (19) 

Thus,  combining  Eqs.  (16)  and  (±7),  CH(x,y;x+£,y+f),  the 
spatial  autocovariance  function  of  the  speckled  irradiance  dis- 
tribution backscattered  from  a rough,  laser-illuminated  target 
observed  through  the  turbulent  atmosphere,  is  given  by: 

>'  CH(x,y;x+e,y+C)  * < H(x,y)  H(x+£,y+C)  > 

!'•;  - < H(x,y)  > < H(x+£,y+£)  > 

( r = < H(x,y)  H(xt£,y+£)  > + < H(x,y)  H(x+4,y+C)  > - « 2 

h „ __  ° _ 

= [h;  + CHs(x,y;x+4,yt0l]  exp  4c^(u, v;x,y;xf£,y+£)  - Hq2 

c s(x,yjx+i,y+0"ir 

1 + — 1 + -2 

H 2 H 2 

o o 
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so  long  as  the  entire  target  lies  within  a single  isoplanatic 
region.  This  equation  may  be  simplified  if  we  assume  further 
that  spatial  stationarity  holds  for  the  separate  autocovariances, 
CH  and  Cjj  (i.  e.,  for  the  speckled  irradiance  in  the  absence 
of  turbulence 2* 3 as  well  as  for  the  turbulence-induced  irradiance 
fluctuations  from  a point  source5*  ®* 10* 12) : 


c~u,o  c„s(e,c)  c„a(4,o  c„s(e,c)  cHau,o 


(20) 


The  spatial  autocovariance  of  the  speckled  irradiance  distribu- 
tion is  thus  the  sum  of  three  terms:  a target-related  laser- 
speckle  term  (Eq.  (17)),  a term  arising  from  the  random  apodiza- 
tion  by  the  turbulence,  and  their  product.  This  result  is  as 
ejected  for  the  correlation  function  of  the  product  of  two 
independent  random  intensity  modulation  effects.18 

Note  that  Eq.  (20)  can  be  solved  in  principle  for  the  auto- 
covariance of  the  target-related  speckle  alone  if  CHa(j-,£)/  Hq2 
can  be  determined  from  independent  measurements. 

Relative  Magnitudes  of  Target-Related  and 
Turbulence-Related  Contributions 


The  first  term  in  Eq.  (20),  CHS(£,£)/  Hq2,  decreases  from 
unity  at  £ = £ = 0 to  zero  when  (£2  + £2)V2  exceeds  a charac- 
teristic "speckle  size".  It  is  clear  from  Eq.  (17)  that  this 
speckle  size  is  inversely  proportional  to  a target-related 
dimension  (e.  g.,  its  maximum  transverse  dimension  if  the  tar- 
get reflects  uniformly  and  is  uniformly  illuminated) . 1-4 
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The  second  term  in  Eq.  (20)  gives  the  covariance  of  the 
turbulence-induced  irradiance  rluctuations  measured  for  a point 
source  with  a pair  of  point  receivers  separated  by  (<;,£).  When 
£ = •£  = 0,  CHa(  (*,£)/  H02  becomes  aH2/  H^2,  which  is  more  com- 
monly written  as  CTj2/  l2  and  called  the  (normalized)  intensity 
variance.10*12  CHa(4,£)/  *Hq2  decreases  from  its  maximum  value 
of  aj2/  I2  at  £ - £ = 0 to  zero  when  (£2  + f^2)1/2  exceeds  a 
correlation  length  of  the  order  of  (AL)1/2,  where  L is  on  the 
order  of  the  path  length  through  uniform  turbulence  (e.  g., 
horizontal  propagation),  10* 19-21  and  on  the  order  of  the  height 
to  the  tropopause  for  vertical  propagation. l0* 19 

For  near -horizontal  propagation  under  "saturation-of-scin- 
tiliation"  conditions, 10 

tfj2/  I2  « 10.  ("Saturation") 

Thus,  except  for  large  values  of  (£,£)  and  speckle  sizes  much 
larger  than  (AL)1/2,  the  covariance  is  dominated  by  the  second 
term  in  Eq.  (20)  under  these  adverse  turbulence  conditions. 

The  observed  irradiance  fluctuations  are  mainly  due  to  the 
scintillation-related  random  apodization. 

On  the  other  hand,  for  vertical  propagation  in  the  visible, 
it  is  known  from  stellar  scintillation  data19  that 

<7j2/  I2  w 0.2.  (Visxble  astronomy;  vertical  path) 

Furthermore,  available  data  at  10.6  njn  wavelength20-22  are  in 
agreement  with  scaling  the  visible  scintillation  data  according 
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to  tT7/"3,  i.  e., 

0j2/  X2  « 0*008.  (10.6  |xm;  vertical  path  at  night23) 

Consequently,  the  autocovariance  of  the  irradiance  distribution 
obtained  at  10.6  jim  from  a near-zenith  exoatmospheric  target 
will  be  dominated  strongly  by  the  first  term  in  Eq.  (20).  The 
observed  irradiance  fluctuations  are  mainly  target-related  laser 
speckles.  The  autocorrelation  function  of  the  speckle  pattern 
is  proportional  to  the  Fourier  transform  of  the  target  auto- 
correlation function,  with  only  a weak  dependence  on  the 
turbulence-induced  random  apodization. 

And  for  near-zenith  targets  in  the  visible  region,  one’s 
ability  to  sort  out  the  contributions  from  each  of  the  three 
terms  in  Eq.  (20)  will  depend  upon  the  relative  sizes  of  (AL)1/2 
and  the  speckle  size,  as  well  as  upon  the  degree  of  turbulence, 
slant  range,  and  zenith  angle.  Based  on  the  measurements  of 
"shadow  bands"  in  stellar  astronomy,  (AL)1/2  can  be  expected  to 
vary  from  3 or  k to  7 or  8 cm  in  the  visible.10  For  small  zenith 
angles,  a j2/  I2  varies  as  sec^G,  where  G is  the  zenith  angle,  and 
p is  a number  on  the  order  of  unity.19 

The  size  of  the  isoplanatic  region  is  difficult  to  esti- 
mate, although  guesses  of  «10~4  rad  for  near-zenith  exoatmospheric 
objects  are  common.  (Even  crude  experimental  data  are  available 
only  for  near-horizontal  propagation. l5* le)  The  present  theory 
thus  is  likely  to  be  appropriate  for  targets  with  transverse 
dimensions  of  meters  at  ranges  of  10  km  or  greater. 
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ition 


By  the  Wiener-Khinchin  theorem, 24  the  power  spectral  den- 
sity of  the  irradiance  distribution  at  the  receiver  aperture 
is  given  by  the  two-dimensional  Fourier  transform  of  the  spatial 
autocorrelation  function  (assumed  to  be  spatially  stationary). 
Thus,  the  power  spectral  density,  S(cd,J2),  is  given  by  the  Fou- 
rier transform  of  <(  H(x,y)  H(x+£,y +£)  )>.  In  the  absence  of 
turbulence,  from  Eqs,  (±6)  and  (17), 

S(cd,£2)  = S = Hq2  6 (to)  6(fl) 

+ I0P  J au  J av  p(u>v)  p(u  • 2T'v  ■ ’ 

—03  —00 

...(21) 


where  6(x)  is  the  Dirac  delta  function,  and 

H^2  is  the  square  of  the  mean  irradiance. 

(The  "lag"  variables,  4 and  £,  are  assumed  to  range  from  -«  to 
This  is  appropriate  when  the  aperture  contains  several 
correlation  "widths" . ) Thus,  in  the  absence  of  turbulence,  the 
power  spectral  density  of  the  received  irradiance  pattern  is  the 
sum  of  a fixed  known  (delta)  function  and  the  target  correlogram. 
i.  e.,  the  spatial  autocorrelation  function  of  the  reflected 
irradiance  distribution  at  the  target. z*  13 

In  the  presence  of  an  intervening  turbulent  atmosphere,  the 
power  spectral  density  may  be  calculated  straightforwardly,  by 
means  of  Eqs.  (l6)-(20),  when  the  entire  target  if  situated  in 
one  isoplanatic  region.  Denoting  results  in  the  absence  of 
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turbulence  by  zero  subscripts,  the  power  spectral  density  is 
given  by: 

s(o>,ft)  - 9 { exp  4c \ * S o(a>,ft)  (22) 

• m 

H 2 ) ° 

o 

where  SQ(ai,fl)  is  given  by  Eq.  (21), 

* denotes  the  convolution  (or  Faltung)  operation,  and 
denotes  the  two-dimensional  Fourier  transform. 

The  convolution  operation  arises  because  the  Fourier  transform 
of  the  product  of  exp  4c^(^,^)  and  the  correlation  function  in 
the  non-turbulent  case  is  given  by  the  convolution  of  their 
Fourier  transforms. 

In  general,  the  power  spectral  density  in  Eq.  (23)  is  the 

sum  of  two  terms.  The  first  term  is  the  uncorrupted  S (cu,  ft), 

o 

which  contains  all  of  the  "fine  structure"  which  may  be  present 
in  the  target  correlogram.  The  second  term  is  a blurred  version 
of  SQ(oD,fi),  obtained  by  convolving  it  with  the  spectrum  of  the 
turbulence-induced  random  apodization.  The  "energy"  contained 
in  this  second  term  is  cHa(0, 0)/  H^2  (i.  e.,  tfj2/  I2)  times 
that  in  the  first  term. 

As  noted  earlier,  for  vertical  propagation  at  10.6  jxm, 

CTj2/  I2  « ±3  so  that  the  first  term  in  Eq.  (23)  is  dominant, 
and  S(cu,ft)  reduces  to  S (o,'ft)  to  a very  good  approximation.  In 
other  cases,  where  turbulent  effects  are,  significant  but  can  be 
determined  from  independent  measurements,  then  the  spectrum  of 
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the  random  apodization  is  known.  Thus,  S can  still  be 

o 

obtained. 

In  any  event,  the  second  term  in  Eq.  (20)  is  determined  by 

the  relative  sizes  of  the  turbulence-related  correlation  length, 

Pq  = (XL)1/2,  and  the  characteristic  speckle  dimensions.  If  the 

uniformly  reflecting  target  is  of  irregular  shape,  the  speckle 

dimension  along  any  given  direction  is  related  to  the  angular 

subtense,  ©t,  of  the  target  along  that  direction  by  dg  w A/©t,  \ 

where  d„  is  the  width  of  the  covariance  function  of  the  ir- 

radiance  distribution,  CHS(£,£),  Eq.  (17).  A feature  related 

to  ©t  will  exist  in  the  target  correlogram,  and  hence  in  the 

power  spectrum,  SQ(oi,f2),  Eq.  (2l),  at  a distance  from  (0, 0)  which 

is  inversely  proportional  to  the  speckle  size,  d Even  in  the 

visible  region,  in  the  limit  that  the  speckle  sizes  are  much 

larger  than  the  scintillation-related  length,  Pq  = (al)1/2, 

s(<o,fi)  in  Eq.  (23)  reduces  to  So(a>,ft),  unchanged.  Physically, 

S (oD,ft)  emerges  in  this  case  because  the  turbulence-induced  ir- 
o 

radiance  fluctuations  are  effectively  aperture-averaged  (i.  e., 
low-pass  filtered)  to  zero  within  a single  speckle  cell.10' 19 
Thus,  target  information  contained  in  the  larger  speckle  sizes 
is  retained  in  the  power  spectral  density,  uncorrupted  by  the 
turbulence . 

In  the  limit  that  the  speckle  sizes  are  much  smaller  than 
p , Eq.  (23)  becomes: 

S(<D,ft)  — —*■  (1  + S-7  I2)  S (o>,fi).  (24) 

Val  » A/et  x 0 
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In  this  limit,  the-  relatively  narrow  spectrum  of  the  scintilla- 
tion imposes  a uniform  apodization  across  the  speckle  power 
spectrum.  The  power  spectral  density  becomes  dependent  upon 
the  variance  of  the  turbulence-induced  irradiance  fluctuations. 
This  variance  itself  is  subject  to  a large  variance. l0> la“22 
Thus,  if  the  speckle  sizes  are  smaller  them  (TJj)1/2,  the 
structure  of  the  power  spectral  density  is  unaffected  by  the 
turbulence,  but  its  overall  strength  is  subject  to  strong 
random  fading  from  look  to  look  in  the  visible.  The  random 
fading  is  negligible  at  10.6  p,m  for  exoatmospherlc  objects  ob- 
served near  the  zenith  at  night,  for  the  factor  multiplying 
So(co,fi)  in  Eq.  (24)  differs  from  unity  by  less  than  1$. 

We  conclude  that,  in  favorable  cases,  so  long  as  the  entire 
target  lies  within  a single  isoplanatic  region,  the  power  spec- 
tral density,  properly  computed,  is  equal  to  the  sum  of  a 
fixed,  known  function  and  the  target  correlogram,  Eq,  (21), 
just  as  in  the  absence  of  turbulence. 
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Summary 

Through  the  use  of  the  extended  Huygens-Fresnel  principle, 
we  have  calculated  the  autocovariance  and  power  spectral  den- 
sity of  the  irradiance  distribution  backscattered  from  a rough, 
laser-illuminated  target  observed  through  the  turbulent  atmo- 
sphere. For  the  special  case  where  the  entire  target  lies 
within  a single  isoplanatic  region,  and  where  the  correlation 
length  of  the  surface  roughness  distribution  is  not  resolved  by 
the  receiving  aperture,  we  verified  that  the  autocovariance  of  the 
irradiance  (normalized  to  the  square  of  the  mean  irradiance)  be- 
comes the  sum  of  three  terms:  the  normalized  covariance  in  the 
absence  of  turbulence  (i.  e.,  the  Fourier  transform  of  the  tar- 
get correlogram) , the  normalized  covariance  of  the  point-source 
irradiance  arising  from  the  turbulence-induced  random  apodiza- 
tion,  and  their  product.  The  power  spectral  density  of  the 
received  irradiance  is  then  the  sum  of  the  uncorrupted,  turbu- 
lence-free spectrum  (proportional  to  the  target  correlogram) 
and  a blurred  version  of  this  spectrum,  obtained  by  convolution 
with  the  spectrum  of  the  random  apodization  produced  by  the 
turbulence.  We  noted  that,  in  principle,  when  the  normalized 
autocovariance  function  of  the  random  apodization  is  known, 
both  the  target  correlogram  and  its  Fourier  transform  may  be 
determined  explicitly  from  the  statistics  of  the  irradiance 
distribution.  Finally,  when  experimental  data  on  stellar  scin- 
tillation were  combined  with  known  data  at  10.6  (xm,  we  concluded 
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that  the  effect  of  the  random  apodization  will  be  negligibly 
small  for  a ground-based  receiver  observing  a suitable  near- 
zenith exoatmospheric  target  illuminated  at  10.6  jim. 
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Appendix  A 

In  this  Appendix,  we  calculate  the  ensemble  mean,  < H(x,y)  )>, 
i,  e.,  the  average  of  H(x,y)  over  all  possible  realizations  of 
the  random  variables,  -t,  /a,  and  /s.  The  ensemble  mean  of  a 
typical  term  in  Eq.  (7)  contains  a factor  of  the  form 

Q(u,v;u»,v*  ;x,y)  = Q(px?pg;p) 

= < exp  [t(px;p)  + -t(p2;p0  cos  H ^(pi“  P£)  *P 

+ ^,s(p1)  - /,s(p2)  + ^(p^p)  - /a(p2;p)II> 

where  p^  = (u,v,0),  pg  = (u‘,v‘,0)  = (u  - 7^,  v - and 

p = (x, y, h),  so  that,  e.  g.,  ^(px_  Pg)*P  = cox  + Qy  . 

If  we  assume  that  -t,  /a,  and  /s  are  effectively  uncorrelated, 
i.  e.,  that  t(pxip)  + ^{p2ip)  is  independent  of  ^ f S(pi)-^* S(P2) 

as  well  as  of  ^(p^jp)  - /a(P2;i>),  then 

Q(px?P2;p)  = < exp  C*’(P1»P)  * -&(i>a»P)U  > Ra  [®XP  ^(Pj/*  P2)*P 

exp  i C^,s(pi)  - ^,S(p2)U>*^  exp  i 0^a(pJL;p)  - /a(P2#pO  ^ 

...(A.l) 

We  shall  show  that,  according  to  our  target  model,7 

< exp  i D**8^)  “ ^'S(P2)U  > = h.  = 6^o  60o  . ...(A.2) 

P1,P2 

Then  (A.l)  becomes 

Q(p.;p  jp)  = Q(p.?p,;p)  6-  _ = < exp  [>(p  jpD  > 

P:.,P2 
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or: 

Q(P1?P2?P)  * < exp  [>t(u,vjx,yO  > 6^  6Qq  ...(A.3) 

Thus,  summing  over  (co,ft)  in  the  ensemble  average  of  Eq.  (7), 
we  obtain  <(  H(x,y)  > = 0,  i.  e.,  Eq.  (8).  It  remains  for  us  to 
show  that  (A. 2)  holds. 

For  f = ft,s(p  ) and  f = ^,s(p  ) real,  jointly  Gaussian 
••  -1  2 2 

random  variables,  it  follows  that 

< exp  i(fi  ± f^)>  = exp  i<  f^>  exp 

• exp  Q-  | Hvar(fi)  + VarCfJ  ± 2 CoV(  f^  f^jj  , 

...(A .4) 

where  Var(fj)  = <(fj  - <fj>^  = < fj2  > - < fj  >2,  j = 1,2, 
cov(fi,f2)  = Cov(f#,fx)  = <(fx  - <fA»(f2  - <f2»>  . 

By  our  assumption7  that  the  correlation  length  between  individual 
scatterers  is  smaller  than  the  resolution  limit  of  the  receiver 

aperture,  we  may  take  Cov(f  , f ) = Var(f  ) _ with  negli- 

12  1 p ^p  • 

gible  error.  Thus,  for  p f p , 1 * 

l a 

< exp  i &'8(px)  ~ ^,S(paO>  = exp  i<  jrf's(px)> 

• exp  -i<  ^,S(p2)>  exp  Q-  ICVarCf^^)  + Var(f2)jJ  . 

s 

Finally,  we  have  assumed  that  the  standard  deviation  of  is 
at  least  equal  to  n, 8 so  that  the  above  expression  is  bounded 

_ 2 _ 

by  e"  n = 5.17  • 10"*5,  effectively  zero.  (A. 2)  follows. 
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In  the  more  general  case  that  p , the  correlation  length 
between  the  individual  scatter era,  is  resolved,  both  the  back- 
scattered  irradiance  and  its  covariance  will  depend  upon 
Cov(f  ,f  J^2*25"2®  Formally,  if  we  assume  a Gaussian  phase 

X 2 

distribution,  so  that  Eq.  ( A-4)  applies,  and  we  suppose  that 

<fx>  83  , Var(fx)  = Var(f2)  = a 2 , and  Cov(fx,f  ) = a 2Cg(Ap*  ), 

, -*  -a  . t \ r r 

where  Ap*  = | p - p | and  C ( Ap* ) is  the  normalized  autocovariance 

12  S 

function  of  the  surface  roughness  distribution, 26  we  obtain  in 
place  of  Eq.  (A.  2) : 

< exp  i Drf*3(p  ) - jrf,s(p  D > = exp  - ["a  2 D-  - C (Ap*  )ZT1  ♦ 

1 2 L s _j 


■ - [»/  & - , 


...(A.5) 


where  C (Ap* ) is  nearly  unity  for  Ap'«  p_,  and  nearly  zero 
s s 

for  Ap*  » p_,  The  additional  complexity  introduced  by  Eq.  (A.5) 

w 

is  unwarranted  in  the  many  practical  cases  in  which  pg  is  un- 
resolved. 2 

However,  it  should  be  pointed  out  that  Var( f . ) may  be  less 

\ A 

than  rr  if  the  target  is  only  slightly  rough  (i.  e.,  imperfectly 
specular).  Then  the  ensemble  mean,  <(  H(x,y)  )>,  may  very  well 
not  be  zero,  and  the  usual  assumption  of  spatial  stationarity2*  25  27 
may  break  down.28 
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By  a "stationary-phase"  argument,  similar  to  that  used  to  obtain 
Eq.  (A. 2)  in  Appendix  A,  we  are  led  to  the  assertion  that,  ac- 
cording to  our  target  model, 7 


< exp  i - jrf,S(jp2)  + ^,S(pg)  - ^'S(P4)H>  ■ 0, 

unless  the  phases  in  square  brackets  cancel  to  zero.  This  can- 
cellation occurs  either  when 

p = p and  p = p (condition  a) 

*1  2 *3*4 

or  when 


p = p and  p - p ( condition  b) . 

1 4 2 3 

For  the  upper  sign  in  (B.3),  we  identify  p = (u,v), 

^ - !“>•  - W"1- 

Then  condition  a corresponds  to  a>  = 0,  = cu'  = ft’  = 0,  and  condi- 
tion b corresponds  to  u = u’,  v = v',  o>  = co’,  '.l  = ft* . For  the 
lower  sign  in  (b.3),  interchange  the  correspondences  of  pQ  and 
p^,  i.  e.,  pg  = (u»,v‘),  P4  = (u‘  - ^j^v*  - Then  the 

condition  a.  again  corresponds  to  u>  = £2  ~ cu*  = £2*  - 0,  and  con- 
dition b corresponds  to  u1  = u - v*  = v - 03  = ~a3'  •> 

0 = -£2» . (These  negative-frequency  possibilities  were  overlooked 
in  Ref.  3,  making  Eq.  (17)  of  that  paper  too  small  by  a factor 
of  2.) 

Substituting  the  non-zero  values  of  Eq.  (B.3)  into  (B.2), 
and  noting  that  Re  exp  (ia)  = Re  exp  (-ia),  we  are  led  immedi- 
ately to  Eq.  (13). 
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As  discussed  in  Appendix  A,  the  additional  complexity  in- 
troduced when  the  correlation  function  of  the  surface  roughness 
distribution  is  taken  into  account  .-explicitly2* 2 ~2  is  not 
warranted  when  the  correlation  length,  p , is  not  resolved  by 
the  aperture, 7 However,  the  calculation  is  straightforward 
if  ^,s(p* ) is  Gaus sian-distributed  and  the  entire  target  lies 
in  a single  isoplanatic  region. 
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FIGURE  CAPTION 

Fig.  i:  An  exoatmospheric  target  is  illuminated  with  coherent 
laser  light.  The  scattered  light  passes  through  the 
turbulent  atmosphere  to  the  ground-based  receiver 
telescope.  The  irradiance  in  the  aperture  is  recorded 
at  each  point  (x,y). 
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Paper  No.  4 


THE  RESIDUAL  EFFECTS  OF  ATMOSPHERIC  TURBULENCE 
ON  A CLASS  OF  HOLOGRAPHIC  IMAGING  AND  CORRELOGRAPHY  SYSTEMS 


(This  paper  was  presented  at  the  Topical  Meeting  cn  Opti- 
cal Propagation  Through  Turbulence.  The  abstract  on  the  first 
four  pages  appears  in  the  Digest  of  Technical  Papers  given  at 
the  meeting. 7 The  Vu-Graphs  used  are  reproduced  here  following 
the  abstract.  Research  Note  N-6/17^-3-33  contained  an  earlier 
version  of  some  of  this  material.3) 

Both  laser  correlography  and  on-axis  lensless  Fourier- 
transform  holography  are  treated  in  this  paper.  The  effects 
of  using  a ground-based  laser  illuminator  are  taken  into  con- 
sideration. C02  laser  correlography  systems  are  shown  to  be 
particularly  insensitive  to  turbulence-caused  degradation  of 
performance  against  space  objects,  although  a valid  verifica- 
tion experiment  cannot  be  performed  over  ground-level  horizon- 
tal propagation  paths. 
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THE  RESIDUAL  EFFECTS  OF  ATMOSPHERIC  TURBULENCE 
ON  A CLASS  OF  HOLOGRAPHIC  IMAGING  AMD  CORRELOGRAPHY  SYSTEMS* 

Michael  Greenebaum 

Riverside  Research  Institute 
New  York,  N.  Y.  10023 

Recent  results  concerning  the  spatial  co- 
variance  of  scintillations  are  used  to  es- 
timate the  effects  of  turbulence  on  holo- 
graphic imaging  and  laser  correlography 
systems. 

In  this  paper,  we  use  recent  experimental  and  theoretical 
results  on  the  spatial  covariance  of  turbulence-induced  scin- 
tillations to  estimate  the  residual  effects  of  turbulence  on 
Fourier-transform  holographic  .imaging  and  on  laser  correlography 
systems.  In  marked  contrast  to  ordinary  imaging,  both  lensless 
Fourier-transform  holography1* 2 and  laser  correlography3* 4 are 
immune  to  the  phase  perturbations  introduced  by  turbulence. 

This  is  true  as  long  as  the  entire  object  lies  well  within  a 
single  isoplanatic  region.1*4  In  particular,  it  is  well  known3 
that  the  two-dimensional  Wiener  spectrum  of  the  speckled  ir radi- 
ance distribution  scattered  from  a coherently  illuminated  object 

■f 
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is  proportional  to  the  ,!correlogram"  of  the  object  illumination 
function.  We  have  shown  previously4  that,  for  a ground-based 
receiver  observing  an  exoatmospheric  object  smaller  than  the 
isoplanatic  patch,  the  phase  perturbations  produced  by  the  tur- 
bulence do  not  affect  this  result.  (The  relevant  isoplanatic 
patch  size  i>3  estimated  to  be  about  10  ^rad  in  angular  extent.) 
The  only  residual  effect  of  the  turbulence  is  a random  apodiza- 
tion  of  the  irradiance.4  Goodman1  has  shown  that,  under  simi- 
lar isoplanatic  patch  conditions,  a similar  result  holds  for 
lensless  Fourier -transform  holographic  imaging.  Thus,  in  each 
of  these  "phase-insensitive"  systems,  estimates  of  the  random 
apodization  are  of  paramount  importance. 

*• 

Because  of  past  difficulties  in  understanding  the  behavior 

of  the  covariance  of  scintillations  in  the  region  of  strong 

turbulence, 5-8  it  has  been  necessary  in  earlier  treatments  to 

restrict  the  system  calculation,  at  least  implicitly,  to  cases 

where  n-2  < 0.3. 4,9  It  has  been  shown  quite  recently10  that 
X ~ 

it  is  possible  to  explain  the  experimental  results  in  the  satu- 
ration region,  > 0.3j  in  a quantitative  manner  by  means  of 
approximate  calculations  based  on  coherence-degradation  arguments. 
In  the  present  paper,  we  use  similar  arguments  to  predict,  at 
least  semi-quantitatively,  the  behavior  of  holographic  imaging 
and  laser  correlography  in  the  region  of  strong  turbulence.  (At 
large  angles  from  the  zenith  direction,  strong  turbulence  will 
be  encountered  even  in  the  infrared.) 
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Secondly,  earlier  analyses1' 2' 4' 9 have  considered  explicit- 
ly the  turbulence  effects  in  the  "down-link"  only,  i.  e.,  from 
the  uniformly  illuminated  exoatmospheric  object  to  the  ground. 
Here,  we  also  consider  the  complications  introduced  when  a ground 
based  laser  illuminator  must  be  used;  i.  e.,  "up-link"  turbulence 
effects  are  included. 

In  particular,  we  can  express  the  resuits  of  our  analysis 
of  laser  correlography  in  terms  of  cus<3^-p)i  the  normalized  spa- 
tial covariance  of  the  irradiances  measured  at  two  points  sepa- 

— *>  # . 

rated  by  distance  p in  the  entrance  pupil  of  the  receiver.  (The 

approximation  of  spatial  stationarity  is  made.)  The  measured 

irradiance,  a stochastic  quantity,  is  related  to  the  transmitted 

laser  power  by  three  random,  independent  amplitude-modulation 

processes,  so  that: 

cusd (?)  - cu(7)  + cs(7)  + ca(p)  + cuQ  c s(7)  + csQ  cs(7) 

+ ca(7)  cu(7)  + cu(7)  cs(T)  ca(7),  ...d) 

where  C (p),  C {~p) , and  C ,(p)  are  the  normalized  covariances  of 
the  received  irradiance  which  would  be  obtained  if  only  one  of 
these  three  processes  were  present  separately.  (The  u,  s,  and  d 
subscripts  relate  to  "up-link"  turbulence,  laser  "speckles",  and 
"down-link"  turbulence,  respectively.)  In  particular,  the  term 
of  direct  interest  in  correlography  is  C (*p),  the  autocovariance 
of  the  "speckles"  which  would  be  observed  in  the  absence  of  tur- 
bulence. If  the  scattering  object  is  a uniformly  reflecting 
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rectangle  of  width  L^,  for  example,  the  characteristic  width  of 
C will  he  Xz/L,  (a  relationship  similar  to  the  van  Cittert- 
Zernike  theorem). 3,4  (The  Fourier  transform  of  Cg  is  the  correlo- 
gram  of  the  object.)  cu(p)  54  cu(o)  is  t^ie  autocovariance  of 
those  irradiance  fluctuations  at  the  receiver  which  arise  from 
fluctuations  in  the  object  illumination  caused  by  scintillation 
and  beam  wander, 6,11  C^(T)  is  the  normalized  autocovariance  of 

turbulence-induced  scintillations  for  a monochromatic  point  source 
at  the  location  of  the  object. 4j 5 3 8j 10  The  behavior  and  influ- 
ence of  and  upon,  one's  ability  to  retrieve  the  object- 

related  information  contained  in  C are  discussed  both  in  the 

s 

weak-  and  strong-turbulence  regimes,  using  techniques  similar 
to  Ref.  10.  It  is  shown  that  a large  receiver  aperture  opera- 
ting with  10-Aim  laser  illumination  can  obtain  relatively  "clean" 
correlograms  of  near-zenith  exoatmo spheric  objects,  and  that 
lensless  Fourier-transform  holography  is  also  feasible. 
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ASSUME: 

1.  OBJECT  LIES  IN  A SINGLE  1SOPLANATIC  REGION. 

2.  "UPLINK'1  AND  "DOWNLINK"  CAN  BETREATED  SEPARATELY. 

3.  OBJECT  CONSISTS  OF  INDEPENDENT,  DIFFUSE  SCATTERERS, 

WITH  RANDOM  PHASE  ANDo*s>X. 

THEN  THE,  INSTANTANEOUS  IRRADIANCE  RECEIVED  IS: 
l(p')  * \$)  i2(p‘)  l3(p') 


m=1  L J 

where  the  three  independent  factors  relate  to: 

TURBULENCE  FROM  TRANSMITTER  TO  OBJECT  (Ij), 

SCATTER  FROM  OBJECT  IN  ABSENCE  OF  TURBULENCE  (l2), 
and  TURBULENCE  FROM  OBJECT  TO  RECEIVER  (1 3). 
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THEN  THE  COVARIANCE  OF  THE  RECEIVED  IRRADI ANCE  IS: 


C,{jo)  = <11  (P‘)  I KP1)!  II (p *+/>)  - 1(P *+P)l> 

l(p')  I (p*  -fyo) 


= Cj^yo)  + C2  (/>)  + Cy(pl 

+ C^(/o)  C2 (/>)  + C2 (/?)  C^lyo)  + c^(p)  C^(yo) 

+ C^fy))  C2 (p)  Cy(p) , 


where 


cm0  ■ 


. <llm(p')  - lm(p'>]  nm(p'+/>)  - lm(p'+/»]> 


1 mtp'  > lm<p'+i°> 


nr 


m = 1,  2,  3. 


As  the  notation  indicates,  we  shall  assume  spatial  stationary, 

i.  e. , the  covariances  Cm  depend  on  />,  not  on  position  in  the 
receiver  aperture,  p'. 
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SCATTER  FROM  OBJECT  IN  ABSENCE  OF  TURBULENCE  (l2): 

Cfp)  = SPECKLE  COVARIANCE 

EFFECTIVE  WIDTH  OF  C20  = Xz/4 

S2(T)  = FOURIER  TRANSFORM  OF  C2£) 

= ' r COR RELO GRAM"  OF  OBJECT 

c2(0)  = 1. 


R- 
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TURBULENCE  FROM  OBJECT  TO  RECEIVER  (13) : 

0,$)  = exp  [4  C/(^)  I - 1 = Cjs(^)/<|>2, 

where  Qf (})  - LOG-AMPLITUDE  COVARIANCE  FOR  POINT  SOURCE 

AT  LOCATION  OF  OBJECT. 


STRENGTH  Oc  RANDOM  APODIZATION 


0. 12  at  0 = 0° 
0.20  at  0 = 40° 


at  0.5/im 


_ 0.003  at  0 r 0° 
0. 005  at  0 = 40° 


at  10. 6yxm  (C02  LASER) 


FOR  ty(0)<  0.3,  NO  PROBLEMS  WITH  "SATURATION". 

Near-zenith  estimates  derived  from  stellar  data  of  Protheroe,  et  al. 
and  scaled  with  wavelength  as  k^: 

C/(0)  ■ 0.028  sec11/6e  (VISIBLE  REGION) 

- 0.0008  sec11/60  (C02  LASER) 

EFFECTIVE  WIDTH  OF  C3£)  « sec  9 (n0  "saturation") 

4 to  10  cm  (VISIBLE  REGION) 

15  to  45  cm  (C02  LASER) 
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SATURATION  REGION >566y i^W 
XXWW'-^ 


' — X - 0.5/xm 

(Protheroe  data  as 
fitted  by  Minott) 


* — X = 10.6/xm 
(Scaling  data  as  X"176) 
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TURBULENCE  FROM  TRANSMITTER  TO  OBJECT  (lx) : 

VARIANCE  OF  OBJECT  ILLUMINATION  IS  LESS  THAN  C3(0) 

BECAUSE  OF  TRANSMITTER  APERTURE-AVERAGING 
(ASSUMING  NO  "SATURATION"). 

GENERALLY  EXPECT  C^)  < C3(0). 

EXPERIMENTAL  VERIFICATION  OF  SCINTILLATION  STRENGTH 
(Minott,  1972). 

C^O)  CAUSES  DYNAMIC  RANGE  PROBLEM  IN  FOURIER-TRANSFORM 
HOLOGRAPHY  WITH  PHOTOGRAPHIC  FILM  DETECTION 
(Goodman,  et  a!. , 1968). 
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NEAR-HORIZONTAL  PROPAGATION: 

"SATURATION"  REGION  ENTERED  AT  LARGE  ZENITH 
ANGLES  WHEN  ty(0)  > 0. 3. 

THEN  1 • C2(0)  « C3(0), 

Cj£>  ACQUIRES  A RESIDUAL  TAIL  EXTENDING  BEYOND^h0  sec  0, 

ANGULAR  EXTENT  OF  ISOPLANATIC  REGION  SHRINKS  WITH 
INCREASING  9.  (Yura,  Clifford  et  ai.) 

EXPERIMENTAL  DATA  FROM  STUDIES  OF  LENSLESS  FOURIER- 
TRANSFORM  HOLOGRAPHY.  (Goodman,  Gaskill,  etal.) 


Example:  For  a 90-meter  horizontal  path,  X=  0.6^.m,  linear 
dimension  of  isoplanatic  region  is  less  than  5 cm. 
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SUMMARY 


CAN  EXTRACT  OBJECT  CORRELOGRAM  (OR  HOLOGRAPHIC  IMAGE) 
FROM  APERTURE-PLANE  IRRADIANCE,  EVEN  WITH  TURBULENCE, 
PROVIDED: 

1.  ANGULAR  SUBTENSE  OF  OBJECT  < 10/xrad. 

2.  OBJECT  ROTATION  IS  NEGLIGIBLE. 

3.  C02  LASER  ILLUMINATION,  DETECTOR  ARRAY,  EVEN  AT  LARGE 

ZENITH  DISTANCES. 

4.  LONG-COHERENCE  VISIBLE  LASER,  OBJECT  VERY  NEAR  ZENITH. 

VERIFICATION  OF  CONCEPT  OVER  NEAR -HORIZONTAL  PATHS  IS 
DIFFICULT. 
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ASSUMPTIONS:  Rough  Object  (crs>\) 

Independent  Scattering  Cells 
Huygens-Fresnel  Principle 
Isoplanatic  Approximation 

INSTANTANEOUS  FIELD  AT  RECEIVER: 

N 

E(p',t)  = 2 E0(p|,t)  G(pj, p* ) , 

H 

where  E0(p,t)cc  r(p)  l0^2(p,t)  exp  iks(p)  , 

G(pj, p* ) - G0(pj,p')  expt^pj.p')  + i^>a(Pj, p')]  , 


l0(p,t)  = illumination  at  wavelength  X= 
r(p)  * local  diffuse  reflection  coefficient, 
sfp)  = random  surface  depth  variable, 

G0(pj,p')  = z*"1  exp  ikz(pj,p')=  Green's  function  with  no  turbulence, 
^(Pj,  p' ) = turbulence-induced  log-amplitude  perturbation, 

<^Pj,p')  = turbulence-induced  phase  perturbation. 


4>a(P|,p')  -<#>a(pj,p')  =^>a(pj,p l+p)  -<£a(Pj,p'+/>)  , 

/{pj,p')  =/(pj,p')  = \ (d,  p'  ) 


(ISOPLANATISM) 


NOTE  IN  F-T  HOLOGRAPHY,  r(pref)  = specular  reflectance, 

s(Pref)  = 0* 
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Paper -No*. ,5. 


COMPUTATIONS  OF  SIGNAL-TO-NOISE  RATIO  IN  CORRELOGRAPHY 


This  paper5  was  originally  written  in  response  to  an  errone 
ous  calculation  of  signal-to-noise  ratio  in  laser  correlography 
which  was  performed  by  the  Hughes  Research  Laboratories.  The 
Hughes  claim  had  been  that  "more  than  10,000  pulses  would  have 
to  be  averaged  to  obtain  adequate  signal-to-noise  ratio"  in 
laser  correlography. 16  In  this  paper,5  it  is  shown  that  only 
20  "looks"  suffice  to  obtain  adequate  signal-to-noise  ratio, 
rather  than  10,000.  Subsequent  laboratory  demonstration  experi- 
ments have  amply  verified  the  signal-to-noise  calculation  given 
here,  as  documented  in  a forthcoming  RRI  Technical  Report. ~5 
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INTRODUCTION 

In  this  note  we  shall  comment  on  the  computation  of  the 
signal-to-noise  ratio  for  correlography  performed  at  Hughes 
Research  Laboratories'1'.  Specifically,  we  show  that  their  dis- 
couraging results  calling  for  10,000  laser  pulses  to  obtain 
adequate  S/N  ratio  at  the  spatial  frequencies  of  interest,  a- 
rises  from  a poor  choice  of  the  processing  used  to  estimate 
the  autocorrelation  function  of  the  speckle  pattern.  (These 
Hughes  predictions  were  obtained  under  the  assumptions  that 
the  speckle  patterns  are  unperturbed  by  the  atmospheric  tur- 
bulence, and  that  the  illumination  energy  is  large  enough  that 
shot  noise  is  negligible-.) 

In  the  following,  we  briefly  review  the  principles  on  which 
correlography  is  based,  emphasizing  the  availability  of  two 
modes  of  processing  of  the  information,  i.e.:  (1)  estimation 

of  the  correlation  function  of  the  speckle  pattern,  and 
(2)  the  estimation  of  the  power  spectrum  of  the  speckle  pat- 
tern. The  statistical,  properties  of  several  estimators  are 
used  to  argue  that  there  exist  statistics  which  can  be  used 
to  extract  intelligence  about  the  object  when  no  more  than 
twenty  laser  pulses  are  used. 
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THE  SPECKLE  .PATTERN 

A coherently  illuminated  object  with  an  optically  rough 
surface  generates  a speckle  pattern.  This  speckle  pattern, 
which  represents  a stochastic  process  over  the  plane  of  the 
receiver,  is  the  "signal"  used  in  correlography.  It  is  well- 
known  that: 

1)  The  autocorrelation  function  of  the  speckle  pattern 
yields  the  absolute  square  amplitude  of  the  Fourier 
transform  of  the  object  ir radiance  distribution. 

2)  The  power  spectral  density  of  the  speckle  pattern 
yields  the  autocorrelation  function  of  the  object 
irradiance  distribution  (correlogram) . 

In  practice,  the  speckle  pattern  I(u,v)  is  available 
only  across  a finite  receiver  aperture,  and  the  number  of 
independent  realizations  of  this  pattern  (obtained  from  a 
series  of  laser  illumination  pulses)  is  limited.  Therefore, 
correlography  is  a technique  in  which  the  information  about 
the  object  is  derived  from  estimators  of  either  the  correlation 
function  or  of  the  power  spectral  density  of  the  speckle  pattern. 
In  correlography,  the  choice  of  the  optimal  processing  method 
must  depend  on  the  statistical  properties  of  these  estimators. 


t 
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CORRELATION  FUNCTION  ESTIMATORS 


The  statistic  used  by  researchers  from  Hughes  to  estimate 

the  correlation  function  of  the  speckle  pattern  is  defined  as: 

1 ^ (n)  _ (n)  _ 

(1)  C^^  (6u,6v)  (I  (u,v)  - I)  (I  (u-6u,  v-6v)  - I)  , 

n=l 

where  M is  the  number  of  independent  realizations  of  the  speckle 
pattern  measured  by  a single  pair  of  detectors  separated  by  the 
distance  (&u,  6v) . The  mean  irradiance, I, is  taken  as  constant 

across  the  pair,  i.  e.,  I = E Cl  (u,v)T=  E Cl(u-6u,  v-6v)] 


The  signal-to-noise  ratio  is  defined  as  the  ratio  of  the 
expected  value  of  C.^  to  the  standard  deviation  of  C^,  and  is 
equal  to 


(2) 


E(Cp 
Var  (C^ 


Vm 


(6u,  6v) 
~=2 


» 


where  C(6u,6v)  is  the  true  value  of  the  correlation  function. 


We  observe  that  if  the  speckle  pattern  has  a negative- 

exponential  probability  density  function  (for  which  the  standard 

deviation  is  equal  to  the  mean)  , C_(^ui_&_y)  ^ for  all  values 

I 

of  (&o,6v). 

* Formula  v-43,  used  to  compute  S/N  ii  the  HRL  report  1 

is  the  square  of  equation  (2)  and  may  be  considered 
equivalent  to  it.  The  reason  for  the  apparent  difference  stems 
from  the  fact  that  v--48  considers  power  ratio  at  the  output  of 
an  electronic  processor,  whereas  our  formulae  (2) , (4)  and  (6) 
express  amplitude  S/N  ratios  for  an  electronic  processor.. 
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Instead  of  the  statistic  C^,  it  is  more  advantageous  to 
introduce  the  statistic  C 2,  which  can  be  realized  by  scanning 
the  irradiance  distribution  across  the  entire  receiver  aperture, 
and  which  is  defined  as 

(3)  C2(6u,6v)  = 7^  * ^ (u,v)  -I^I^  (u-6u,v-6v) -I^jdudv, 


where  A is  the  area  of  the  receiver.  The  signal-to-noise 
ratio  is  defined  as  before,  i.e.,  as  a ratio  of  the  expected 
value  of  C2  to  its  standard  deviation.  In  the  particular  case 
of  a circular  Gaussian  autocorrelation  of  the  speckle  irradiance, 
for  example,  signal-to-noise  is  given  by 


(4) 


Speckle  size 


Speckle  size 


C (6u,6  v) 
i2 


t 


where  it  is  assumed  that  the  true  value  of  the  correlation 
function,  C,  is  given  by 


0 = I2  e-2h2Ao2  _ 
where  h^=  6u^+6v^ 

and  ho/  vr  is  the  1/e  width  of  the  correlation  function. 

Comparing  the  statistics  and  C2  > we  observe  that  al- 
though both  are  unbiased  estimators  of  the  correlation  function  C, 
only  the  second  one  is  a consistent  estimator.  Thus,  as  the 
aperture  of  the  receiver.  A,  increases,  C2  tends  to  be  clustered 
about  the  ensemble  average. 
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The  signal-to-noise  ratio  computed  for 

is  superior  by  a factor  / A ^ with  respect  to  that 

Y Speckle  area 

computed  for  C^.  For  a 5-m-diameter  target  at  a distance  of 
400  km  from  the  ground-based  receiver,  illuminated  with  0-5  pm 
laser  light,  and  observed  with  a receiver  of  1 -m  aperture  diameter 

this  factor  is  on  the  order  of  25 ! (or  625  when  the  square  of  the 
Eq4  is  compared  with  formula  v-48) . 

POWER  SPECTRAL  DENSITY  ESTIMATORS 

2 3 

Based  on  rough  analytical  predictions  ’ , and  as  sup- 
ported by  laboratory  demonstrations  at  RRI4,  we  have  found 
that  the  power  spectral  estimators  of  the  speckle  pattern  have 
even  more  useful  properties,  and  are  likely  to  be  the  optimal 
statistics  for  oorrelography. 


One  such  statistic,  which  has  been  used  successfully  to 
obtain  an  adequate  signal-to-noise  ratio  with  less  than  twenty 
independent  samples  of  the  speckle  pattern,  is  defined  as 


(5)  Sl(f.  *.,  -i  • | jf 


1^  (u,  v)  e 


i2n(ftlu+fv  v) 


dudv 


The  signal-to-noise  ratio  computed  for  this  statistic  is 
equal  to 
(6) 


S1  (VarS1)  2 


IN 
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at  any  spatial  frequency  "resolvable"  by  the  receiver  aperture 

A.  In  deriving  the  last  equation,  it  has  been  assumed  that  the 

2 3 

estimator  S^,  has  a chi-square  distribution  ' with  two  degrees 
of  freedom  (for  which  the  mean  is  equal  to  the  standard  devia- 
tion) . For  M = 20,  the  signal-to-noise  is  equal  to  4.5  for  all 
arguments  of  the  power  spectrum  estimator. 
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Paper  No.  6 


PRELIMINARY  PLAN  FOR  A LASER  CORRELOGRAPKY  EXPERIMENT 


(This  paper  is  a slightly  edited  version  of  RRI  Research 
Note  N-3/174-3-35* 2 It  describes  the  design  of  a feasibility 
demonstration  experiment,  using  coherent  ruby-laser  illumination 
of  suitable  space  objects,  and  designed  around  a more-or-less 
conventional  telescope  facility.  It  is  in  the  form  of  Vu-Graphs, 
with  commentary  on  facing  pages,  and  expands  upon  some  of  the 
issues  discussed  in  Sections  II  A and  B of  Paper  No.  1.) 
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The  purposes  of  this  experiment  are  threefold.  First 
and  primary  is  the  demonstration  of  the  feasibility  of  this 
technique  for  obtaining  correlograms  of  orbiting  space 
vehicles.  Second,  a meaningful  opportunity  to  assess  the 
practical  utility  of  such  signatures  will  arise  from 
successful  completion  of  the  experiment.  Finally,  several 
aspects  of  laser  correlography  system  design  may  be  confirmed 
by  the  experimental  results.  These  include  the  influences  of 
glint,  target  motion,  atmospheric  apodization  and  choice  of 
optimum  signal  processing. 
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This  figure  shows  the  basic  elements  of  a laser 
correlography  system  configuration.  A ground-based  laser 
illuminates  a space  vehicle;  the  backscattered  laser  light 
forms  a self-interference  pattern  which  is  sampled  by  the 
irradiated  receiver  array.  It  is  expected  that  the  structural 
nature*  of  this  pattern,  which  is  related  to  the  shape, 
size  and  structure  of  the  target,  will  survive  passage 
through  the  relatively  thin  turbulent  atmospheric  layer. 


* More  precisely,  the  two-dimensional  power  spectrum  of  the 
pattern.  See  Paper  No.  3 for  the  gory  details,  which  demonstrate 
that  this  power  spectrum  is  proportional  to  the  correlogram  of  the 
target. 
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A laser  correlogram  is  the  autocorrelation  of  the  image 
of  the  object.  Such  signatures  should  be  useful  for  measure- 
ment of  various  object  dimensions  as  well  as  the  identification 
of  features  on  targets  and  entire  target  systems.  Resolution 
is  expected  to  be  limited  almost  totally  by  receiving  aperture 
size  with  minor  degradation  in  the  visible  region  where 
the  proposed  feasibility  experiment  will  be  performed. 

Finally,  it  should  be  noted  that  several  looks  or  returns 
are  required  from  a single  target  aspect  in  order  to  obtain 
measures  of  the  average  speckle  structure;  our  restate  have 

to  be  adequate  for  reducing  the  variability 
(See  Paper  No.  5.) 


shown  thirty  looks 
of  the  signature. 
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AVERAGE  SPECKLE  STRUCTURE 
30  FRAMES  ADEQUATE 
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This  block  diagram  shows  the  essential  elements  that 
must  be  specified  in  order  to  describe  the  laser  correlography 
experiment. 
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BLOCK  DIAGRAM  OF  LASER 
CORRFLQGRAPHY  SYSTEM 


CORRELOGRAM 
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For  the  purpose  of  experimental  design,  the  target  was 
assumed  to  be  described  by  the  parameters  of  the  accompanying 
figure.  Also,  atmospheric  seeing  of  better  than  2 arc  seconds 
and  transmission  loss  of  1.5  dB  were  assumed. 
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Seeing  ~ 2 

Atmospheric  two-way  loss  = 1.5  dB  at  0.7/i.m 
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This  graph  plots  the  average  return  as  a function  of 
target  range  for  10-  and  25-joule  laser  pulses  directed  in 
5 arc  second  beam  widths  at  targets  described  above.  The 
left-hand  ordinate  shows  the  average  number  of  photons 
arriving  per  detection  cell  per  transmit  pulse;  the  right- 
hand  ordinate  reduces  this  number  according  to  the  6%  quantum 
efficiency  of  an  S-25  photocathode.  In  each  case,  the  area 
of  a detection  cell  is  equal  to  one -fourth  the  area  of  the 
smallest  speckle  in  the  backscatter  pattern.  Note  that  a 
10- joule  pulse  directed  at  a target  at  350-km  range  provides 
8.5  photoelectrons  per  detection  cell.  It  is  expected  that 
the  actual  laser  energy  per  pulse  will  lie  somewhere  between 
the  10  and  25  joule  curves  shown. 
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As  a result  of  this  calculation,  the  laser  illuminator 
specification  expressed  in  tne  accompanying  chart  is  obtained. 
A nominal  10- joule/pulse  ruby  transmitter  is  specified, 
and  the  required  coherence  length  of  10  meters  allows  for 
5-meter  target  depth. 
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Also,  as  a result  of  the  above  calculation,  the 
receiver  intensity  detector  must  be  signal  shot  noise  limited 
and  our  studies  show  that  this  can  be  achieved  with  a 
dynamic  range  of  25  dB  using  a SIT  or  EBS  tube  preceded  by  an 

image  intensifier. 
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Finally,  the  output  resulting  from,  e.  g. , a two- 
minute  satellite  pass  will  be  30  to  50  frames  of  data  with 
target  resolution  of  20  cm  and  a signal-to-noise  r&? io  of 
~3.  During  a two-minute  pass,  it  is  necessary  that  the 
target  have  a nearly  constant  range  and  small  rotation;  such 
targets  are  available. 
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